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I. INTRODUCTION 


Combustor liner durability is one of the major challenges in high pressure ratio turbopropulsion engines. 
To resolve conflicting combustor performance requirements, along with enhanced structural durability in 
a cost-effective manner, conventional (empirical) combustor design techniques need to be complemented 
by multidimensional aerothermai analysis. 

To improve predictive capabilities of aerothermai models, NASA Lewis Research Center sponsored 
Phase I aerothermai modeling activities (Srinivasan et al, 1983; Kenworthy et al, 1983; Sturgess, 

1983*) to assess current state-of-the-art numerical schemes and physical models. The main objectives of 
aerothermai modeling Phase II are: 

1 . Develop advanced numerical schemes 

2. Collect benchmark quality experimental data to quantify interaction between dome swirlers 
and primary jets 

3. Collect benchmark quality data for fuel nozzles and their interaction with dome swirlers 

4. Use advanced numerics and benchmark quality data to validate advanced aerothermai models 

Two advanced numerical schemes were developed by Karki et al (1988) under NASA HOST (Hot-Sec- 
tion Technology) sponsorship. The experimental activities and model validation efforts for nozzle- 
swirler interaction were conducted jointly at University of California at Irvine and Allison Gas Turbine 
Division of General Motors Corporation (Nikjooy et al, 1993). This report summarizes a joint Allison/ 
Purdue University effort on benchmark quality data and model validation for swirler-jet interaction for 
an idealized nonreacting primary zone. 

This report consists of seven sections, including Section I, and three appendices. 

In Section II, an overview of gas turbine combustor flowfield characteristics is presented. This is fol- 
lowed by a description of an idealized combustor flow model and an explanation of the integrated mod- 
eling/ experimental approach. 

Section III presents the details of the experimental rig and instrumentation. It begins with a discussion 
of initial rig design and operation, and followed by a description of rig refinement activities and final 
test rig. Finally, flow visualization, instrumentation, and data reduction results are presented. 

In Section IV, the experimental data are presented and discussed. This includes dome annular jets only, 
primary jets only, and dome annular jets and primary jets. 

In Section V, the problem of calculating turbulent flows is posed more precisely by introducing and dis- 
cussing the averaged equations governing the mean flow quantities. The appearance of turbulent trans- 
port terms in these equations makes apparent the necessity of introducing turbulence models. The mod- 
els are discussed in order of increasing complexity. The details of the solution procedure adopted for 
the highly coupled and nonlinear governing equations are explained next. 

In Section VI, the computational results are comparted with experimental data. The results are pre- 
sented for the standard k-e model, algebraic second-moment (ASM) closure, and differential second- 
moment (DSM) closure. 

Finally, Section VII summarizes the main conclusions that emerged from this study and puts forward 
some recommendations for future work. 


* References for Section 1 are listed at the end of the section. 
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II. SELECTION OF EXPERIMENTAL CONFIGURATION 


2.1 GAS TURBINE COMBUSTOR FLOW FIELD CHARACTERISTICS 

Gas turbine combustion systems need to be designed and developed to meet many mutually conflicting 
design requirements, including high combustion efficiency over a wide operating envelope and low NOx 
emissions, low smoke and low lean flame stability limits and good starting characteristics, low combus- 
tion system pressure loss, low pattern factor, and sufficient cooling air to maintain low wall tempera- 
ture levels and gradients commensurate with structural durability. The flow field around and within 
the combustor liner (Figure 2.1-1* ) is quite complex in that it includes swirl, regions of recirculation, 
fuel injection, atomization, fuel evaporation, mixing, turbulent combustion, soot formation/ oxidation, 
and convective and radiative heat transfer processes. The phenomenological understanding of these 
processes is not well established, and the relevant nonlinear coupled transport equations are difficult to 
solve. 

The combustor design and development process has been empirically based with limited help from mul- 
tidimensional calculations. A number of correlations have been used by combustor designers to help dur- 
ing the design and development activities. Many researchers (Lefebvre, 1984 and 1985; Plee and Mel- 
lor, 1978 and 1979 + ) have proposed semiempirical correlations for gaseous emissions, smoke, lean 
blowout, ignition, pattern factor, and combustion efficiency. Professors Lefebvre and Mellor have de- 
veloped very useful correlations that can be used for scaleup, for data correlations, and for providing 
some insight. These correlations or their variants are being used by the gas turbine industry (Steele et 
al, 1987; Rizk and Mongia, 1989). 

The empirical/analytical combustor design methodology introduced by Mongia and Smith (1978) has 
been used for designing a number of gas turbine combustors (Mongia, 1982; Mongia et al, 1986). The mul- 
tidimensional calculations provide a good understanding of combustor internal flow field and therefore 
can be used for guiding a combustor design process. However, because of incomplete understanding of the 
various combustion processes and numerical diffusion, the three-dimensional calculations for practical 
gas turbine combustors cannot be considered quantitatively accurate (Srinivasan et al, 1983; Kenworthv 
et al, 1983; Sturgess, 1983). Some progress is being made in further improving the numerical methods 
(Karki and Mongia, 1989), however more effort is needed to achieve the capability required for accu- 
rately predicting combustor performance parameters including radial profile, combustion efficiency, 
smoke and gaseous emissions, and wall temperature levels and gradients. 

Although encouraging qualitative comparisons between data and calculations were achieved for the 
practical combustors (Rizk and Mongia, 1991), it became quite dear that a significantly increased level 
of effort is required to achieve the model accuracy required for providing definitive guidance during de- 
sign process. Moreover, each important process of combustion (e.g., turbulence, kinetics, turbulence/ 
kinetic interaction, spray, etc) should be investigated separately and in combination so as to improve 
fundamental understanding. Model calculations were, therefore, performed for existing data from 
simple to complex flow under HOST sponsorship (Srinivasan et al, 1983) and the following main 
conclusions were made: 

1. Improve phenomenological understanding of nonreacting flows relevant to those encountered in 
gas turbine combustors. 

2. Make benchmark quality measurements in critical areas of interest. 

3. Undertake a systematic model validation effort to identify areas of further improvement. 


Figures for Section II appear at the end of each subsection. The figure number identifies the subsection 
in which the figure is discussed. 

+ References for Section II are listed at the end of the section. 
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Figure 2.1-1. Model predictions of flow field around and within combustor. 
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2.2 IDEALIZED COMBUSTOR FLOW MODEL 


The overall combustor flow field can be divided into two categories, jet flows and swirling flows. Each 
of these subflow categories has been extensively studied but, surprisingly, little attention has been fo- 
cused on the combined problem. 

2.2.1 let Injection into a Cross Stream 

Jet injection into a cross stream has been widely studied. A paper by Crabb et al (1981) includes a list of 
24 earlier experiments. Many of the early studies dealt with single jets flowing into unbounded streams, 
so the results were not directly relevant to primary or dilution jets in combustors. In real configurations, 
the jets interact with each other and the blockage causes the main stream to accelerate. Walker and 
Kors (1973) and Holdeman et al (1973) have reported measurements of the temperature distribution 
downstream of a row of cold jets injected from one side wall of a rectangular duct into a hot cross flow. 
This configuration was more relevant to gas turbine combustor design than previous studies, although 
the jet-to-mainstream density ratio was somewhat lower than what prevails in practice. Holdeman 
and Walker (1977) and Cox (1975) analyzed these experimental results and derived correlations that 
can be used for preliminary design of dilution jets in annular combustors. Although these correlations 
are used in practice, the failure of the original experiments to accurately simulate the flow conditions 
of a combustion chamber reduces the accuracy of the predictions. For example, the cross flow was essen- 
tially isothermal whereas the flow at the plane of the dilution jets in a practical combustor has both 
radial and circumferential temperature gradients. Second, the cross flow was not swirling, so these ex- 
periments more closely simulate the conditions that prevail for dilution jets than for primary jets. The 
swirl and recirculation created at the dome strongly influence the penetration and mixing of the pri- 
mary jets. Finally, the fine details of the jet injection hole can influence the mixing process, and this 
boundary condition was not accurately simulated. Since no velocity measurements were made and the 
upstream boundary conditions were not determined, assumptions of numerical solutions must be consid- 
ered when these results are compared with experimental results. 

Similar experiments explicitly intended to model the injection and mixing of dilution jets have also 
been performed (Srinivasan et al, 1982). In addition to studying the effect of jet spacing and diameter 
and the jet-to-mainstream density and momentum-flux ratios, the experiments also included more real- 
istic effects such as nonuniform free-stream temperature distribution and convergence of the duct. 

The results for isothermal flow in a constant area duct agreed very well with the earlier data of 
Walker and Kors. The jet spreading characteristics in the nonisothermal cases were very similar to the 
isothermal case, but the temperature distribution was, not surprisingly, found to be strongly determined 
by the original upstream temperature distribution. For isothermal flow in a convergent duct, the jet 
mixing rates in both the radial and transverse directions were generally greater than in the constant 
area duct. However, when the free stream temperature was made nonuniform, the details of the distri- 
bution strongly influenced the mixing. When the higher temperature was toward the wall with the 
jets, the lateral and transverse spreading rates were not appreciably altered from the isothermal re- 
sults. When the higher temperature was along the wall opposite the jets, the spreading of the jets was 
significantly slower. 

To analyze the data, the correlations derived by Holderman and Walker (1977) were simply modified 
to account for nonuniform temperature distributions and area change. The new results gave good quali- 
tative predictions, but needed refinement to be quantitatively useful. These data are useful for prelim- 
inary design of the dilution holes, but do not apply to the primary zone, largely for the same reasons 
discussed in connection with Walker s experiments. 

Srinivasan et al (1983) did an extensive review of measurements of simple and complex turbulent flows 
that could be used as benchmarks for numerical simulations of the aerothermodynamic phenomena that 
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occur in gas turbine combustion chambers. Computations of many of these flows using advanced codes, 
and with several turbulent transport models, were also carried out. 

Kamotani and Greber (1972) have measured velocity and temperature distributions in a heated jet in- 
jected into a semi-infinite cross stream. Mean velocity and turbulence intensity were measured with a 
constant current anemometer, a yaw probe was used to determine the direction of the mean velocity vec- 
tor, and total temperatures were measured with a small thermocouple. The trajectory of the jet was de- 
scribed in terms of the loci of the points of maximum velocity and temperature. The jet momentum ratio 
was the major factor that controlled the trajectories, although the temperature centerline was weakly 
influenced by the density ratio. The path of peak velocity was always located farther from the wall 
than the line of peak temperature. Two counterrotating streamwise vortices were generated on either 
side of the jet centerline which strongly influenced the subsequent behavior of the jets. 

Crabb et al (1981) reported detailed turbulence measurements in the same flow. Crabb's jet was the 
same temperature as the main stream, but helium was injected in the jet to act as a tracer for jet fluid. 
Measurements of the mean and root mean square (rms) fluctuation velocities in the streamwise and nor- 
mal directions were made near the jets with a single component laser Doppler velocimeter (LVD) sys- 
tem. Downstream, where the turbulence intensity was less than 0.30, hot wires were used to measure all 
three velocity components (mean and rms) and two Reynolds shear stresses. Probablility distributions 
were determined with both instruments. A gas sampling probe was used to determine the mean concen- 
tration of helium. The measurements showed that away from the injection plane, the locus of maximum 
velocity does not correspond to jet fluid, even though the initial jet velocity was larger than the free 
stream velocity. This explained Kamotani's results. The Reynolds stress distributions indicate that 
gradient diffusion models for momentum, such as the k-e model, will not yield accurate results, particu- 
larly in the wake region near the wall. While this single-jet flow is only a subcomponent of the jet- 
swirl interaction, both of these sets of measurements could be used for some benchmark computations. 

Rho et al (1988) measured the turbulent characteristics of a cross jet mixing after the geometric cross 
point of the centerlines of the two nozzles. The flow has been generated by two circular nozzles and the 
crossing angle of the jets was fixed at 45 deg with respect to each other. The flow field is strongly influ- 
encedby the cross point stagnation region, which creates an elliptical jet region. The measurements 
have been performed at Reynolds numbers of 5.2 x 10 4 and 6.5 x 10 4 and some influences of nonsimilarity 
have been studied. 

Catalano et al (1989) reported measurements and computations for the flow of a turbulent jet discharg- 
ing into a cross flow confined between two parallel plates. For jet-to-cross flow velocity ratios equal to 2 
and 4, mean and fluctuating velocity components were measured by laser Doppler anemometry. Results 
show that impingement of the jet on the opposite wall occurs at a velocity ratio of 4. 

Koutmos and McGuirk (1989) reported measurements of mean velocity components and the corresponding 
turbulence intensities in a water model of a can-type gas turbine combustion chamber. The effect of 
variation of the flow split between the swirler and the dilution holes on the flowfield in the primary 
zone was investigated. At low levels of swirler flow rate, a large toroidal recirculation with high lev- 
els of turbulence energy were generated in the case of the combustor. As the swirl level increases, the 
strength of this recirculation was observed to weaken. The flow pattern downstream of the secondary 
jets was the same for all flow splits due to the strong mixing. 

Sherif and Pletcher (1989) reported measurements of the velocity and turbulence characteristics of a 
round turbulent isothermal jet in cross flow at velocity ratios of 2, 4, and 6. Their experiments were ear- 
ned out in a water channel 8.53 m long and 0.61 m x 1.067 m in cross section. Mean velocities, Reynolds 
stresses, and correlation coefficients were measured using anemometry methods. 
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Liou et al (1989) reported studies on the three-dimensional jet-jet impingement flow in a closed-end 
cylindrical duct with two 60 deg side inlets. The impingement of jets provided an effective way to mix 
the fluids. The measurements were made using the laser-Doppler velocimeter. The Reynolds number of 
the air flow was 2.6 x 10 4 . Simple geometrical expressions were deduced and were found to effectively 
estimate the stagnation points and the optimal head heights. 

Vranos et al (1991) conducted an experimental investigation of cross stream injection and mixing. The 
transverse injection from slanted slot and orifice injectors were studied. The main parameters of their 
study were momentum-ratio, density ratio, and injector geometry. Measurements of injectant concentra- 
tion distributions were obtained. The data indicate that a density ratio greater than unity retards mix- 
ing. Furthermore, above a certain momentum-flux ratio, mixing with slanted slot injectors was better 
than with round hole injectors. 

2.2.2 Confined Swirling Flows 

Interest in confined swirling flows is very strong because of the highly intense turbulence that can be 
created. If the swirl is sufficiently strong, a recirculation region is created which can act as a flame- 
holder. Most gas turbine combustion systems use swirl to establish a primary combustion zone in the 
dome region, where the equivalence ratio is close to one. Syred and Beer (1974) and Lilley (1977) re- 
viewed the state of knowledge of swirling reacting flows almost a decade ago. Experimental and ana- 
lytical advances have since increased the body of knowledge for both reacting and nonreacting flows. A 
review of experimental investigations of swirling confined flows was also done by Srinivasan et al 
(1983) as part of a larger assessment of the accuracy of the aerothermodynamic submodels. Both simple 
and complex flows were considered and computations using advanced numerical codes and several turbu- 
lence models were carried out. Published measurements by Escudier et al (1980), Altgeld et al (1983), 
Rhode et al (1983), Vu and Couldin (1980), and Janjua et al (1982) were reviewed. The most important 
conclusions drawn were that the k-e model and the corresponding gradient diffusion scalar transport 
model worked well in simple flows but were not satisfactory for recirculating swirling flow. Algebraic 
stress models yielded similar results from the mean velocity and scalar fields but better predictions of 
the Reynolds normal stresses. The Reynolds shear stress predictions were not satisfactory. Similar con- 
clusions were drawn by Habib and Whitelaw (1980) in an experimental and analytical study of con- 
fined swirling coaxial jets. 

Rhode et al (1983) measured mean velocity fields in a cylindrical chamber supplied from one end by an 
annular swirling jet. Chambers with two inlet geometries, a sudden expansion and gradual (45 deg) ex- 
pansion, were studied. Attention was focused on the effects of increasing amounts of swirl on the recircu- 
lation regions and on the mean velocity profiles. Numerical predictions of the flow were presented in a 
companion paper (Rhode et al, 1982) and Janjua et al (1982) measured the turbulent characteristics of 
swirling and nonswirling flows in the same geometry. One of the interesting observations was that 
swirl causes a recirculation zone near the center line in which axial velocities are very low but the cir- 
cumferential velocity is maintained. 

As previously mentioned, these particular measurements were used by Srinivasan et al (1983) to evalu- 
ate the accuracy of current numerical codes. However, because no jets were injected into the flow, these 
measurements have limited direct application to gas turbine combustor design. A recent paper by Fer- 
rell et al (1984) considers the case of injection of radial jets into a cylindrical duct. Both numerical and 
experimental results were presented. The k-e model was used to predict a flow of six jets into the duct. 
Swirling and nonswiriing flows were computed and it was found that radial penetration of the jets is 
reduced by the swirl. The swirl was sufficiently weak that no recirculation was induced. An experi- 
mental investigation of the flow of a single radial jet into a nonswiriing flow was also described. Flow 
visualization photographs and mean velocity profiles only were presented. The experiments focused on 
the effect of jet velocity ratio and a downstream area reduction on the jet penetration. Experiments such 
as these are beginning to approach the complexity of the primary zone flow. Further continuation of 
this work to the measurements of the turbulent quantities will provide valuable results for the combus- 
tion designers. 
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Gouldin and his coworkers at Cornell University have published a series of investigations of reacting 
and nonreacting swirling flows. Vu and Gouldin (1980) used pitot and hot-wire probes to investigate the 
nonreacting flow in a cylindrical chamber fed from one end by coaxial jets. The strength and direction of 
the swirl in each jet could be independently controlled. In all of the cases reported, the average axial 
velocity of the outer jet was about 0.7 times that of the inner. Both coswirling and counterswirling flows 
were studied, and the strength of recirculation zone was found to increase as the swirl in the outer jet 
was reduced from positive (coswirling) to negative (counterswirling) values. This was associated with 
increased turbulent dissipation within the interjet shear layer. Emissions measurements made in a ge- 
ometrically similar combustion chamber were directly correlated with these results (Oven et al, 1979). 
The counterswirling flow had lower combustion efficiency because the vigorous mixing in this shear 
layer rapidly quenched the combustion reactions. 

In a recent paper, Halthore and Gouldin (1984) reported measurements of mean and fluctuating density 
made in the recirculating primary zone of a cylindrical combustion chamber. Methane at an equiva- 
lence ratio of 0.80 was supplied through a central swirling jet located on the axis of the chamber. More 
air was supplied through an outer coaxial jet that could also be made to swirl. Time-resolved measure- 
ments of total molecular number-density were made using Rayleight scattering and the mean number- 
density of methane molecules was measured with Raman scattering. The results showed that the recir- 
culation zones created by counter and corotating jets are different. Velocity measurements in both 
nonreacting and reacting flows in the same combustion chamber have also been reported (Gouldin et al, 
1983). Together these data represent another flow that can usefully serve for verification of numerical 
codes, but again, because no primary zone jets were injected, the data have limited direct applicability 
to the design of gas turbine combustors. 

Johnson and Bennett (1984) and Roback and Johnson (1983) have made very careful and complete mea- 
surements of confined coaxial jets. Their geometry was basically the same as that of Habib and 
Whitelaw (1980) and they too studied both swirling and nonswirling flows. The most important aspect 
of this work is that data on both turbulent mass and momentum transport were taken. Velocity mea- 
surements were made with one and two component LDV systems. The concentration of fluorescent dye 
initially injected only in the central jet was determined by measuring the intensity of the fluorescence 
from a small volume of fluid illuminated with a laser. Under appropriate conditions the intensity of 
the fluorescence is directly proportional to the number-density of fluorescent molecules. By observing 
the crossing volume of a dual beam LDV system with two photodetectors, one of which used a filter 
which passed only the laser radiation while the second used a filter which passed only the fluorescent 
radiation, turbulent transport of the dye, defined by the cross correlation between velocity and concen- 
tration fluctuations, was measured. The velocity concentration data pairs were analyzed to determine 
probability density functions for axial, radial, and circumferential velocity fluctuations, concentration 
fluctuations, and their various products. Moments of these distributions up to the fourth-order were cal- 
culated. In addition, detailed surveys of the velocity distributions at the exits of the swirlers were 
measured. This very complete data set represents an excellent example of a benchmark experiment 
which can be used for validating current and advanced prediction codes. Again, however, the lack of 
secondary jets reduces the direct applicability of the data to the understanding and design of flow in 
the dome region of a combustor. 

Many excellent studies of various aspects of the fluid dynamic phenomena in gas turbine combustors 
have been performed. Most of these, however, have omitted key factors which would have greatly in- 
creased the applicability of the results throughout the combustor. Jet injection studies have omitted 
swirl and recirculation in the primary stream, thus restricting their usefulness to dilution jets near the 
combustor exit. Studies of swirling recirculating flows have typically omitted jet injection, but most gas 
turbine combustion systems require additional combustion air to be injected through holes in the dome to 
establish the correct fuel/air ratios for stable combustion. 
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23 INTEGRATED MODELING/EXPERIMENTAL APPROACH 


An experimental rig of rectangular cross section was selected for this investigation (Figure 23-1). The 
rig satisfies two important requirements, it simulates many important features of advanced annular 
combustors and it has been designed so that all of the boundary conditions will be accurately known. 

The annular combustor model was designed to simulate the cold flow characteristics of current gas tur- 
bine combustors. Air enters the test section through five swirlers mounted in the head plate at one end 
of the rectangular duct. The swirler center-to-center spacing is equal to the duct height, and the vanes 
are flat plates set at an angle of 60 deg with respect to the centerline of the nozzle. Air through two 
separately controlled lines enters through the primary jets located in the top and bottom of the duct. 
The 3.0 in. x 15.0 in. (73 cm x 38.1 cm) cross section has the same hydraulic diameter as a typical ad- 
vanced combustor. 

A 3-D code was used to predict the flow field for a number of test configurations to help identify the ex- 
perimental test matrix. The calculations determined areas of strong velocity gradients where mea- 
surements would be taken. Since the purpose of the flow interaction program is to obtain benchmark 
data with which to verify 3-D turbulence models, it is necessary to involve the 3-D model in the selec- 
tion of the experimental configuration. This is important for several reasons. First, the 3-D model can 
help identify significant features of the configuration, and in so doing can help determine which pa- 
rameters are important and which are not, thus concentrating effort on fruitful test procedures. Sec- 
ondly, the 3-D model is also useful in establishing boundary conditions and the overall scope of the 
work. 

Before the proposed experimental configuration was selected, a number of preliminary cases were com- 
puted using a 3-D code. The calculations were made for the two-equation k-e model. There were sev- 
eral reasons for not using advanced turbulence models. First, the models are only illustrations used to 
help select the experimental configurations. Secondly, the results were obtained with a relatively 
coarse grid. Since these calculations are not grid independent, there is an excessive and unknown amount 
of numerical diffusion, thus obscuring the advantages offered by advanced models. 

Five preliminary cases were run with a grid size of 22 x 14 x 14. This was considered sufficient to de- 
scribe each case adequately without requiring excessive computer time. Using a finer grid would have 
increased the accuracy but would not have affected the trend of the results. 

The first case (designated as baseline) consisted of parallel axial swirlers side-by-side in a rectangular 
duct (in a manner analogous to annular combustor design) with opposing primary jets in line with the 
swirlers and downstream one duct height (3.0 in.). The swirler spacing was 3.0 in, and the swirlers had 
a tip diameter of 1.5 in. and a hub diameter of 1.0 in. The primary jet flow was set equal to 75% of the 
swirler flow, the flow rates were scaled for inlet velocities of 400 ft/sec. 

Velocity vector plots at four longitudinal planes are shown in Figure 2.3-2. The planes are located at 
0.5 in. increments, starting midway between swirlers with the last plane, in line with the center of the 
swirler. These plots show the flow patterns as they are affected by the swirler and the primary jets. 

Similar velocity vector plots at six cross-sectional planes are shown in Figure 23-3. These are located 
at distances of 0.75, 1.5, 235, 3.0, 3.75, and 4.5 in. (x/H = 0.25, 0.5, 0.75, 1, 1.25, 1.5) from the swirler, re- 
spectively. From these plots it can be seen that the swirl quickly degenerates into two transverse 
streams until the streams are broken up by the primary jets. One interesting result is that the transverse 
flow direction is reversed by the primary jets. 

Figure 2.3-4 shows contour plots of the axial velocity and the kinetic energy of turbulence in the longi- 
tudinal plane through the center of the swirler. The axial velocity contour plot shows the shape of the 
recirculation zone extending from the swirler hub to the primary jets. The neck in the recirculation zone 
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suggests that the primary jets should be located nearer to the swirler to obtain a more efficient recircu- 
lation flow. The contours of kinetic energy of turbulence also suggests the same conclusion. 

In the second case, the primary jets were moved closer to the swirler to a distance of one-half test sec- 
tion height, i.e., 1.5 in. Analogous plots are shown in Figures 2.3-5, 23-6, and 23-7. The contour plots of 
Figure 2.3-7 indicate that the recirculation zones set up by the swirler and by the primary jets have coa- 
lesced and the overall recirculation is reduced. A primary jet location somewhere between 1.5 to 3.0 in. 
would probably be better. 

The third case was similar to the baseline except that the primary jet flow was doubled. The plots for 
this case are presented in Figures 23-8, 2.3-9, and 23-10. Again, the axial velocity contour plots of Fig- 
ures 2.3-9 and 2.3-10 indicate some coalescing of the recirculating region. For primary jets at this loca- 
tion the baseline flow should be increased a little less than double from the baseline. 

The fourth case was an annular version of the baseline configuration. For this case the inner wall radius 
was 5.5 in. and the outer wall radius was 8.5 in. The analogous plots (Figures 2.3-11, 23-12, and 2.3-13) 
show a deflection of the flow toward the outer wall. While there are differences between the annular 
and rectangular configurations, the similarities of the trends make the rectangular case a worthwhile 
test of the computer model. 

The last case was the same as the baseline but without the primary jets. The plots for this case (Figures 
2.3-14, 2.3-15, and 2.3-16), when compared with those for the baseline case, indicate that the effect of 
the primary jets is quite apparent. 
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Figure 2.3-1. Nominal distribution of survey points within the test section. 
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Figure 2.3-2. Longitudinal velocity vector plots for baseline configuration. 
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l-igiire 2.3-3. Cross-sectional velocity vector plots for baseline configurations. 
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Figure 2.3-4, Contour plots for plane through center of swirler for baseline configuration. 
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Figure 2.3-5. Longitudinal velocity vector plots for configuration with primary jets at x/H= 0.5. 
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Figure 2.3-6. Cross-sectional velocily veclor plots for configuration with primary jets at x/H= 0.5. 
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Figure 2.3-7. Contour plots for plane through center of swirler for configuration with primary jets of x/h 

= 0.5. 
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Figure 2.3-8. Longitudinal velocity vector plots for configuration with doubled primary jet flow. 
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Figure 2.3-9. Cross-sectional velocity vector plots for configuration with doubled primary jet flow. 
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Figure 2.3-10. Contour plots for plane through center of swirler for configuration with doubled primary 

jet flow. 
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Figure 2.3-1 1. Longitudinal velocity vector plots for annular configuration 
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Figure 2.3-12. Cross-sectional velocity vector plots for annular configuration. 




Figure 2.3-13. Contour plots for plane through center of swirler for annular configuration. 
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Figure 2.3-14 Longitudinal velocity vector plots for configuration with swirler only. 
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Figure 2.3-15. Cross-sectional velocity vector plots for configuration with swirler only. 
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Figure 2.3-16. Contour plots for plane through center of swirler for configuration with swirler only. 
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HI. EXPERIMENTAL RIG AND INSTRUMENTATION 


The experimental phase of the work was performed at Purdue University. An experimental rig that 
was considered adequate for obtaining "benchmark” quality data was designed and constructed by Alli- 
son Gas Turbine for use by Purdue University. It was anticipated from the outset that the jet-swirl in- 
teractions would be extremely complex and possibly unsteady, with many aperiodic features. Accord- 
ingly, it was decided that an experimental rig that could be operated with water as the working fluid 
and that would be suitable for exploration of various features of interaction through flow visualization 
would be a valuable adjunct facility to the air-operated facility to be utilized for obtaining quantita- 
tive data. Such a rig, identical in all respects to the air-operated rig, was designed and constructed by 
Allison Gas Turbine and provided to Purdue University for investigations. 

Preliminary investigations with the air-operated and the water-operated rigs confirmed the existence 
of various dynamical features in the fiowfield. It also became clear that the experimental results 
would be extremely sensitive to relatively small and apparently insignificant differences in the geo- 
metrical features and initial flow conditions in the two test rigs. After a substantial amount of testing, 
and analysis of fiowfield data and observations, it was concluded that at least the following three fea- 
tures needed careful investigations: (i) the geometrical veracity of the test section; (ii) the symmetry 
of the swirler that was utilized to generate the swirling flow; and (iii) the fiowfield generated by the 
primary jets, including both velocity distribution and unsteadiness of flow. A series of small-scale ex- 
perimental studies was undertaken on each of the aforementioned flow features by building individual, 
dedicated experimental rigs and carrying out both quantitative and flow-visualization tests. The re- 
sults obtained with such rigs are of considerable interest in themselves. They revealed both the type of 
complex interactions arising in the fiowfield and the accuracy required in the construction of the two 
large air- and water-operated rigs. 

Based on such data, the original rigs were modified for the conduct of experiments and the generation of 
the main body of fiowfield data. 

It should be pointed out that the preliminary studies conducted on the two test rigs also revealed the 
difficulties in utilizing standard intrusive and nonintrusive measurement techniques and means in such 
complex flowfields. Thus, it was found that the use of laser Doppler velocimetry (LDV) in the air-op- 
erated rig, and that of nephelometry in the water-operated rig, required considerable development in 
application. The final body of test data provided was obtained following all of the desired modifica- 
tions to the basic test rigs and to the means of diagnostics and measurement. 

The data presented in this section are extracted from M.S. theses written by Barron, 1986; Morgan, 1988; 
Seal, 1988; and Dwenger, 1990, respectively. 

3.1 INITIAL RIGS - DESIGN AND OPERATION 

The results of extensive experimentation in the initial air and water rigs are presented in the next two 
subsections. 

3.1.1 Five-Swirler Air Rig 

The annular combustor model was designed to allow measurements in a fiowfield similar to that exist- 
ing in current annular combustors. These measurements are to be used as a data base in the development 
and analysis of computational methods for the design of future gas turbine combustors. 

The annular combustion chamber model, several views of which can be seen in Figures 3.1. 1-1 and 
3.1. 1-2*, was designed to simulate the cold (nonreacting) flow characteristics of a current gas turbine 
combustor; the characteristics were swirl, recirculation, primary cross jets, and high levels of turbu- 

Figures for Section III appear at the end of each subsection. The figure number identifies the 
subsection in which the figure is discussed. 
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lence. The recirculation zone, typical of mixing zones in current combustors, is generated by the swirl 
imparted to the inlet annular jet flow by a set of swirl vanes. The length of this recirculation zone is de- 
termined by the location and strength of the primary cross jets. In order to allow the best possible mod- 
eling of actual combustor conditions, the annular and primary cross jets were fully developed prior to en- 
try into the test section. 

In an ideal annular combustion chamber geometric similarity exists for each swirler. This geometric 
similarity is in the form of a volume, or cell, of the flowfield extending downstream of the swirler inlet 
plane. The flow patterns and characteristics in each of these cells will be the same, along with the 
boundary conditions. The arrangement of these cells can be seen in Figure 3.1. 1-3. Therefore, in order to 
understand the characteristics of the flowfield there needs only to be consideration of the region or ceil 
corresponding to one swirler. 

The same principle is utilized in the experimental apparatus. Due to the possibility of cell-to-cell 
mass transfer resulting from secondary flow patterns, as occurs in actual combustors, the decision was 
made to use multiple cells in the experimental apparatus. By allowing the adjacent cells to interact 
the development of these secondary flows was possible. The choice of five cells over three cells was 
made on the basis of the two extra cells adding more interaction between cells and, hence, more realism. 

In an actual combustor a cell would be slightly curved, as can be seen in Figure 3. 1.1-4. Comparing this to 
the combustor model cell, it can be assumed that the effect of the curvature is slight or negligible, thus 
allowing the duct and the swirler cell to be constructed with a simple rectangular cross section. 

The duct rectangle cross section has a aspect ratio of 5, with a duct height of 3.00 in., a width of 15.00 in. 
and a cross section corresponding to each swirler cell of 3.00 in. by 3.00 in. The duct has a downstream 
extension from the swirler headpiate of 10 duct heights, or 30.00 in. For this experiment the configura- 
tion consisted of two opposing primary cross jets mounted along the cell centerlines. The centerlines of 
the jets were located one duct height, 3.00 in., downstream of the test section inlet plane. 

The duct was originally constructed with 0.500 in. thick plexiglass plates throughout. This choice of 
materials proved unsuitable for some of the LDV measurements due to factors such as: 

• loss of beam crossover at the probe volume 

• loss of pinhole alignment 

• change in effective focal distance to the receiving lens 

All three of these factors can be attributed to the slight variations in the thickness of the plexiglass 
sheets. The first factor can also be attributed to scratches in the surfaces of these sheets. Scratches are 
a problem with plexiglass applications such as this due to the need for continuous cleaning of the sur- 
faces. 

From a series of preliminary diagnostic studies, it was found that the measurements that were the most 
dramatically affected by the use of plexiglass were those made in the xy plane. These measurements 
were affected more than those in the xz or yz planes due to the fact that the beams passed through the 
duct sidewalls at a considerably larger distance from the probe volume. For the measurements in the xz 
or vz planes this distance ranged from 0.40 in. to 2.60 in., measured from the inner wall, while the that 
for the xy plane measurements ranged from 6.00 in. to 9.20 in. Since the beams must travel farther, the 
effective change in the beam angle due to thickness variations or scratches becomes more pronounced, re- 
sulting in an increase in the signal loss. 

The solution to this problem was the replacement of the plexiglass sidewalls with 0.25 in. thick glass 
plates. These plates facilitated unimpeded beam traversal due to very uniform thickness as well as 
high resistance to scratching. 
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The annular jets were constructed out of stainless steel, with a length of 9.55 in. and inner and outer 
radii of 1.084 in. and 1.459 in., resulting in a length-to-effective diameter ratio of 25. The swirlers were 
actual engine hardware currently being used in the Allison 570-K turboshaft gas turbine engine (Figure 
3. 1.1-5). These swirlers, also constructed from stainless steel, are vane-type with 12 flat vanes at blade 
angles of 60 deg. The approximate swirl number is calculated to be 1-52. A diagram of the annular jet 
can be seen in Figure 3.1. 1-6, while a separate diagram of the swirler can be seen in Figure 3.1. 1-7. 

In order for the experimental apparatus to be representative of an annular combustor, the five annular 
jets must pass the same mass flow. Since the pressure drop across a jet can be related to the mass flow, or 
average velocity, through the jet 


Ap = ypuL g 

where 


( 1 ) 


U = mass averaged velocity in the jet 
k = loss coefficient of the jet 
Ap = pressure drop across the jet 
Li^vg = m / (Ap) 

Two pressure taps were installed on each jet, the distance between the taps being as large as geometri- 
cally possible to allow the measurement of the largest pressure drop. The larger the measured pressure 
drop the smaller the percent uncertainty of the calculated velocities. 

Control of the pressure drop, and the corresponding mass flow, is allowed by the throttle valves 
mounted at the inlets of the annular jets. These throttles, which also can be seen in Figure 3. 1.1-6, are 
infinitely adjustable to allow either an increase or decrease in mass flow. 

Initially, the swirlers were mounted flush, as in Figure 3. 1.1-8, but preliminary flow visualization 
showed that upon exiting the swirler the flow attached to the headplate in the manner shown. Figure 
3. 1.1-9 shows the modified installation of the annular jets and swirlers to the headplate of the duct. 
The swirlers are recessed a small distance from being flush with the headplate. 

This type of flow phenomena is not desirable in gas turbine combustors since the fluid would transfer a 
large amount of heat to the combustor liner with the main flow losing energy. The increased heat trans- 
fer rate, as well as increased surface shear stress due to the larger velocity gradients near the wall, 
would serve to reduce the life of the combustor, increasing the maintenance costs. The loss in energy 
from the main flow results in reduced available turbine inlet temperature and a corresponding loss in 
power. The desired type of swirler exit flow pattern can be seen in Figure 3.1. 1-9. This pattern produces 
a highly stable flame and regions of high combustion intensity. From these flow visualization results, 
the swirlers were recessed in the annular jets a distance of 0.25 in. in order to produce the desired flow 
pattern at the test section inlet. 

The preliminary flow visualization was performed in the single swirler apparatus which, as stated 
previously, operated at a smaller annular pressure drop than the model combustor. Further flow visu- 
alization in the single swirler apparatus of the test section inlet flow pattern, with the swirler in- 
stalled in the recessed position, was also performed. These studies indicated that a bistable flowfieid 
pattern existed for a small annular jet pressure drop range very near the maximum of the apparatus. 

For pressure drops (AP) of less than 0.0178 psi, the desired flow pattern existed, while for pressure 
drops of greater than 0.0190 psi the attached flow pattern existed. In the range 0.0178 psi < AP < 0.0190 
psi, the flow pattern was bistable, switching between the two established flow patterns in a seemingly 
random fashion. From these flow visualization results it can be assumed that the model combustor, op- 
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erating at a pressure drop of 0.040 psi, exhibited the attached flow pattern. Due to time constraints, 
further modifications of the test section inlet to produce the desired flow pattern were not possible. All 
the measurements were, therefore, conducted with the attached inlet flow pattern. 

The primary cross flow jets were constructed from stainless steel with a length of 10.75 in. and an inside 
diameter of 0.43 in., also resulting in a length-to-diameter ratio of 25. As with the annular jets, the 
primary jet mass flow balance was maintained by balancing the pressure drops on the jets. The pressure 
taps were mounted to give the largest possible pressure drop and the inlets of the jets were fitted with 
throttle valves to allow the tailoring of the pressure drops. A diagram of a sample primary cross jet 
and throttle valve can be seen in Figure 3.1.1-10. 

An important characteristic in the analysis of a combustor of this type is the primary jet-to-annular jet 
mass flow ratio. Originally, the desire was to have a primary jet mass averaged exit velocity of 300 
ft/sec and an annular jet mass averaged exit velocity of 150 ft/sec. These velocities gave the desired 
mass flow ratio of 0.81, a total duct mass flow of 0.517 Ibm/sec and a duct average velocity of 21.6 ft/ sec. 

Although the five annular and ten primary cross jets appeared identical to each other, some inconsis- 
tencies did exist. Differences in the maximum measured velocities were discovered during preliminary 
testing. Since the jets were required to have the same mass flow, the maximum pressure drop was lim- 
ited to the lowest measured maximum of all the jets. 

The maximum average velocity for the annular jets was 93.4 ft/sec, while for the cross jets 250.0 ft/sec 
was the maximum. These measurements were made with a pitot tube and the LDV. The original mass 
flow ratio could have been obtained by adjusting the throttles valves on the cross jets to reduce their 
mass flow; the decision was made to have the highest possible jet exit velocities in order to keep the 
turbulence levels as high as possible. The resulting mass flow ratio was 1.08, with a total duct mass 
flow of 0.371 lbm/ sec and a duct average velocity of 15.8 ft/ sec. 

In order to draw the ambient laboratory air into and through the test section some type of fan or blower 
had to be utilized. A 30-hp Chicago Blowers centrifugal blower was selected, due to its ability to pass 
the large required mass flow. Connecting the test section duct to the fan was a large volume plenum 
chamber. The large volume was chosen to minimize the possible effects of feedback of larger scale fluc- 
tuations into the test section. Preliminary tests indicated that some large-scale feedback in the form of 
large-scale velocity fluctuations did exist, prompting the use of screens and straws at the end of the test 
section. The additions did not cause a significant loss in mass flow and did eliminate the previously ob- 
served feedback. The plenum chamber and the centrifugal fan can be seen in Figure 3.1.1-11. 

3. 1.1.1 Inlet Conditions 


Flow profiles showing the development of the primary jet flow into the test section can be found in Fig- 
ures 3.1.1-12 through 3.1.1-27 (note, y = 1.5 in. is the horizontal centerplane, z = 7.5 in. is the vertical 
centerplane, and z = 6 in. and z = 9 in. are the vertical midplanes). These profiles, showing only the 
vertical (y direction) velocity component (corresponding to the jet axial velocity), are measured at the 
vertical locations y=0.06, 0.25, 0.50, 0.75, 1.00, 1.25, and 1.50 in. These mean and root mean square (rms) 
velocity distnbutions should show profiles similar to pipe flow with substantial variations with 
vertical distance. 

The measurement grids in both the xy (same as yz) and xz planes can be seen in Figures 3.1.1-28 and 
3.1.1-29. The xy plane grid shows a change in the number of grid points in the x direction, with a corre- 
sponding change in the z direction, at the third vertical station from 11 to 13 points, which was insti- 
tuted since the jet flux would spread at the larger y distances. The xz plane grids corresponding to these 
various y locations also show the change in the grids with Figure 3.1.1-29 (a) illustrating the grid for 
y=0.06 in. and y=0.25 in. and Figure 3.1.1-29 (b) illustrating the grid at y=0.50, 0.75, 1.00, 1.25, and 1.50 
in. Both of these grids use a spacing of 0.05 in. in both the x and z directions. 
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At the y=0.06 in. station a close comparison to classic pipe flow data such as that from Laufer (1954*) 
would be expected. A comparison of Laufer’s data for both the mean and rms quantities, at Reynolds 
numbers (Re) (based on maximum velocity and jet diameter) of 50,000 and 500,000, to the data from the 
current experiment (Re = 75,000) can be seen in Figures 3.1.1-14 and 3.1.1-15. The data from both Laufer 
and the combustor model are nondimensionallzed by the centerline (maximum) velocity. From Figure 
3.1.1-14 it can be seen that the experimental V profile from the combustor model shows better agreement 
with the Re=500,000 data than the Re=50,000 data. Slight differences are to be expected, such as the 
slight spread of the flow due to the dissipation of the sharp gradients near the jet wall, since the com- 
bustor model measurements were made 0.06 in. from the jet exit. 

Figure 3.1.1-15 shows plots of turbulence intensity (the rms velocity divided by the centerline mean ve- 
locity) for both the combustor model data and Laufer's two cases. As can be seen, both of Laufer’s cases 
show strong peaks near the walls of the jet with maximum values near 9%. The combustor model data 
show a maximum near 13%, but the peak could be larger due to the rather poor resolution of the mea- 
surement grid (compared to that of Laufer). The centerline value for the combustor model is also larger, 
by nearly a factor of two, than both of Laufer's data sets. This large difference between the data sets is 
mainly attributed to the effect of the opposing cross jet as well as the effect of the local main flow. 

As discussed previously, flow visualization of the jet stagnation point region indicated that the loca- 
tion of the stagnation point fluctuated, in a rather random fashion, within a 0.75 in. range between the 
exits of the two jets. The region between the two jets is evidently highly turbulent and unsteady, result- 
ing in high levels of turbulence intensity across the test section, even near the cross jet exits. Since the 
centerline of the cross jets was located at z=7.50 in., the jet flow issued into a region of cross flow of ve- 
locity of roughly 10% of the maximum jet velocity. Although this is a relatively small magnitude 
crossflow velocity, this effect also resulted in an increase of the turbulence intensity. 

At the y=0.25 in. station (Figure 3.1.1-16) the V distribution shows a slight spread in the jet flux, with 
not much reduction in the maximum velocity. The Vnns velocity shows a large increase in the turbu- 
lence fluctuations over those found at y=0.06 in. The values measured at this location in the combustor 
model are much larger than the values that would be expected from pure jet flow, again due to the op- 
posing cross jets creating a region with very high levels of turbulence. 

At y=0.50 in., the V distribution can be seen to be translated in both the x and z directions, this skewness 
corresponds to a bending of the centerline of the jet flow due to the mam flow acting as a cross flow and 
exerting a shear force on the jet flux. The spread of the jet flux is increased over at y=0.25 in., with the 
maximum velocity remaining near 325 ft/sec. The distribution shows overall increases across the 
grid, with the peak levels, relatively sharp at the y=0.25 in. location, being dissipated to give a 
smoother contour. 

The V distributions at y=0.75, 1.00, and 1.25 in. indicate a decrease in the maximum vertical velocity, 
as well as continued spread in the vertical jet flux, with distance from the jet exit. The location of the 
center of the vertical jet flux can also be seen to show little translation in the xz plane. The Unus dis- 
tributions show the decay of the circular peak region into a concentrated peak at the y=1.25 in. station 
with an xz plane location at the center of the vertical jet flux. This peak occurs in the stagnation point 
location fluctuation region, resulting in a very high level of turbulence. 

The y=1.50 in. plane is located in the middle of the stagnation point fluctuation region and, thus, shows 
high levels of turbulence as at y=1.25 in. with the rms velocity contour displaying a much smoother dis- 
tribution than the previous stations. The main point of interest is the reduction in the vertical jet flux, 
the main fluid motion in this region is in the xz plane as seen in the xy and yz plane vector plots. Also, 
from these vector plots the location of the stagnation point was illustrated to be near y=1.35 in.; from 


References for Section III are listed at the end of the section. 
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these results the stagnation point is not evident even at y=1.50 in. Since the different data sets for the 
primary jet and for the main flow were measured at different times, the operating conditions could 
have been slightly altered. The peak vertical velocity has been reduced from 180 ft/sec, at y=1.25 in., 
to 50 ft/sec, indicating that the stagnation point is close to the y=1.50 in. location. 

3. 1.1.2 Mean Flowfield 

Three-dimensional information on the main flowfield has been obtained through the use of a single- 
component LDV system. Distributions of U, V, W, U', V', W', LTV', and K are calculated from the indi- 
vidual data samples. 

The choice of a measurement grid depends on the character of the fluid motion. More detailed mea- 
surements are needed in regions of large gradients and complex flow patterns, with less detail necessary 
in regions of simpler, less complex flow. Due to the three-dimensionality of the flowfield, a choice was 
made as to the regions of the flow in which to concentrate the measurements. The yz plane was in need 
of a measurement grid with greater spatial resolution than that for the xy plane, due to the smaller 
scale structures to be found in the yz plane. Thus, measurements were to be made in yz planes for various 
downstream locations, with these downstream locations chosen so as to also be more concentrated in the 
more complex regions. 

Another consideration in measurement grid sizing involves the probe volume size, the resolution of the 
grid should be such that probe volume overlap from one data point to the next does not occur. Rotation 
of the test section was necessary to allow the measurement of the W characteristics. Due to this rota- 
tion, the length of the probe volume, on the average being 0.05 in., is the driving factor for the grid 
spacing in both the y and z directions which was chosen to be 0.20 in., or 4 probe volume lengths, for both 
directions. 

The spacing in the x direction of the various planes was chosen as 0.50 in. in the region of the more com- 
plex flow, 0.0 in. < x < 4.00 in., with larger spacing thereafter. The downstream locations of the yz 
measurement pianes were x=1.00, 1.50, 2.00, 2.50, 3.00, 3250, 4.00, 4.50, 6.00, and 9.00 in. The yz plane and 
xy plane measurement grids can be seen in Figures 3.1.1-30 and 3.1.1-31. 

Although a complete scan in the yz plane from one wall to the other was desired, limits on the vertical 
travel of the probe volume were imposed. The range of vertical probe volume travel was, therefore, 
limited to 0.40 in. < y < 2.60 in. Since data pertaining to only one cell was desired, the limits on the 
travel of the probe volume in the z direction were 5.80 in. < z < 9.00 in., where the z=5.80 in. location 
allows comparison to a geometrically similar station at z=8.80 in. 

In the measurement of the W characteristics some geometrical restrictions were encountered at x=2.50, 
3.00, and 3.50 in., necessitating the removal of the cross jet throttles and mounting plate, very little re- 
sulting flowfield alteration occurred. The measurements at the x=2.5Q in. and x=3.50 in. locations were 
then able to be made while those at x=3.00 in. were not, due to the necessary removal of the cross jet 
supporting structure. 

The mean flowfield found in the model combustor is best represented by vector plots in the xy and yz 
planes. These vector plots allow the relative values of the velocities at points in the flowfield to be 
seen thus illustrating the general mean flow characteristics. The xy plane vector plots can be found in 
Figures 3.1.1-32 through 3.1.1-37 while the yz plane vector plots can be found in Figures 3.1.1-38 through 
3.1.1-43. 

From the xy plane vector plots two mam characteristics can be seen, the formation of a general toroidal 
recirculation zone and the influence of the primary cross jets. A great amount of interaction exists be- 
tween the recirculation zone and the cross jets. As a consequence of this interaction, the recirculation 
zone and the cross jet interaction must be discussed together, not separately. 
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The recirculation zone is generated by the application of swirl to the inlet annular jet flow and the step 
change in area at the inlet to the test section. This zone is created in order to enhance the combustion 
process. The general character of this zone can be seen in the xy plane vector plots. This region of back- 
flow exists over a large portion of the cell cross section, showing asymmetric behavior. 

The purpose of the primary jets is to shorten the recirculation zone in order to decrease the overall 
length of the combustor. This goal has been achieved as no backflow exists past the x=3.00 in. station, 
but while the length of the recirculation region is shortened by the addition of the jets, the strength of 
recirculation has been enhanced by the radial jet effect from the cross jet stagnation point. The location 
of the stagnation point is controlled by the relative strength of the two opposing jets, if the jets had 
equal mass flows (pressure drops), the stagnation point would be located in the vertical center of the 
duct. From Figures 3.1.1-34 (a) and 3.1.1-35 (a), it can be seen that this location is not in the vertical 
center of the duct, but at roughly y=1.35 in. Preliminary flow visualization had illustrated large-scale 
fluctuations of the location of the cross jet stagnation point in the range 1.00 in. < y < 2.00 in., indicating 
the unsteadiness present in the flowfield. 

The influence of the cross jet stagnation point on the recirculation zone, and the downstream flow, is eas- 
ily seen in Figures 3.1.1-34 and 3.1.1-35. The flow can be imagined to flow out of the jets, reach the stag- 
nation point, and be forced to move in a radial-tvpe motion in the xz plane. The recirculating backflow 
thus receives a driving force which generates stronger recirculation. 

Previous studies on similar experimental apparatuses (Brondum and Bennett, 1986; Rosenweig et al, 
1961; McLaughlin and Tiederman, 1973) show a recirculation region which is symmetric about the duct 
centerline in the xy plane. These studies, though, used axisymmetric (cylindrical) test section geome- 
tries, whereas the geometry used in this study was of rectangular cross section, with a square cell cross 
section. As movement is made equal distances away from the centerline location, z=7.50 in., with a 
square cross section expectations would be for the vector profiles to be the inverse of one another. Since 
all of the cells are assumed geometrically similar, a necessary condition is that the flow conditions at 
the plane boundaries separating the cells be the same. At the z=6.00 in. and z=9.00 in. boundaries, the 
flow conditions should be the same. 

Figure 3.1.1-38 illustrates the regions of recirculating flow in the yz plane at x=1.00, 1.50, 2.00, and 2.50 
in. The recirculation zone at x=1.00 in. appears to be located off-center in the positive z direction, with 
strong recirculation existing into the upper right comer and a smaller amount of recirculating mass flow 
exists in the upper left comer. This backflow existing in the upper left comer could be part of the recir- 
culation zone from the adjacent cell. From the vector plots for 7.40 in. < z < 8.00 in. the effect of the at- 
tached inlet flow pattern is evident, the recirculation is strong across the duct height with positive ax- 
ial flux occurring only near the walls. Strong positive axial flux can be found to exist in the lower left 
and lower right comers. The profiles at the z boundaries are roughly similar in form, with larger mag- 
nitudes existing at z=9.00 in. 

The vector plots indicate the structure of the flow patterns at x=1.50 in. to be very similar to that at 
x=1.00 in. The recirculation can be seen to extend the entire width of the cell cross section with a larger 
portion of this backflow existing in the upper y range. Nearer the z=6.00 in. and z=9.00 in. boundaries, 
the lower y range is dominated by strong positive axial flux. The profiles at the z boundaries, like 
those at x=1.00 in., are similar in form but with much larger magnitudes existing at z=9.00 in. 

Except in the range 7.00 in. < z < 8.40 in. the vector profiles at x=2.00 in. and x=2-50 in. are very similar. 
Near the center of the ceil, the strong influence of the cross jets on the recirculated mass flow is illus- 
trated. The extent of the backflow in the cell cross section can be seen to be reduced for the x=2.50 in. lo- 
cation than for the x=2.00 in. location, the backflow is more concentrated near the center of the cell cross 
section. This can be attributed to the flow from the cross jet stagnation point region. The strongest influ- 
ence of the stagnation point on the recirculated mass flow is in the center of the cell cross section. For 
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both the x=2.00 in. and x=250 in. stations, the vector plots on the z boundaries are very different in 
shape and magnitude. 

The behavior of the mean flowfield for stations downstream of the cross jet injection plane is relatively 
simple compared to that of the upstream stations. The main characteristic is the influence on the local 
flowfield from the radial jet flow pattern formed at the cross jet stagnation point. This effect is best 
seen in Figures 3.1.1-34 and 3.1.1-35. As one moves in the z direction, either positive or negative, the 
vertical spread of the radial jet can be seen. This spread is an effect of viscous dissipation and is strong 
enough that, for the x=3.50 in. plane, no influence of the cross jets can be seen at either the z=6.00 in. or 
the z=9.00 in. boundary. The dissipation of this effect occurs in the axial direction, as well, since no 
vertical variation of any consequence in the axial velocity profile can be seen at x=6.00 in. or x=9.00 in. 

The yz plane vector plots (Figures 3.1.1-38 through 3.1.1-43) allow the illustration of the flow devel- 
oping in the various cross-stream planes. At all of these stations some type of smaller scale swirling 
motion exists, with the largest concentrated near the cross jet injection location. Unexpected, however, 
was the absence of overall swirl. Some amount of overall swirl was expected in the first duct height of 
the test section due to the large amount of swirl (swirl number: S=1.52) applied to the inlet flow. 

As with the vector plots in the xy plane, the vector plots in the yz plane show the effect of the cross jets 
on the flowfield. The far downstream locations, x=6.00 in. and x=9.00 in., show very little motion in the 
cross-stream plane, illustrating the large amount of dissipation present. 

The vector plots of most interest are those of the x=2.50 in. and x=3.50 in. downstream locations. At both 
of these stations, the strong effect of the cross jets on the combustor model flowfield can be seen. The 
flow pattern in the yz plane can be seen to be dominated by four cells of vortex-type motion. For both of 
these locations, the left side of the plane (looking downstream, z < 7.50 in.) shows a much more sym- 
metric (with the y=1.50 in. centerline) pattern than the right side. 

The driving force behind this type of motion is, as in the case of the xy plane, the flow resulting from 
the cross jet stagnation point region. The flow can be seen exiting the cross jets, reaching the stagnation 
region, and flowing in a radial type motion in the xz plane. This radial motion exhibits itself in the yz 
plane as fluid motion in the positive and negative z directions. Thus, with generally no cross flow 
between the adjacent swirler cells, this jet effect drives this vortex-type cell motion. 

Farther upstream of the x=2.50 in. station one can still see the swirling motion, although not nearly as 
structured as the x=2.50 in. or x=3.50 in. cases. The x=2.00 in. location shows two regions of vortex-type 
motion, with a more structured flow pattern possibly altered by the existence of the large amount of 
cross-flow at the z=9.00 in. boundary. This type of cross-flow is not seen at the z=6.00 in. boundary, and 
thus is a probable cause of the asymmetry in the flowfield. This cross-flow is unexpected since all the 
swirler cells were assumed to be geometrically similar, resulting in similar flowfields for the different 
cells. This cross-flow is partially responsible for the asymmetries found in the right half of the grid at 
the x=2.50 in. and x=3.50 in. locations. 

Cross-flow can be seen to exist at the x=1.50 in. plane and, to a lesser extent, at the x=1.00 in. plane. 
These two locations share similarities in flow patterns other than aforementioned cross-flow from the 
adjacent cell. In each of the vector plots of these two locations, two regions of vortex-type motion can be 
observed, one in the upper right corner, and a second, much weaker, vortex motion located near the mid z 
and lower y ranges. The behavior in the lower z range for the full y range at both the x=1.00 in. and 
x=1.50 in. stations is rather undefined. At x=1.00 in., scattered or bad data seem to obscure flow details, 
while at the x=1.50 in. station, no prevalent type of motion exists. 

As stated earlier, no overall swirling motion is found in the first duct height of the test section, as 
would be expected from the large amount of swirl applied to the inlet flow. The combination of the 
primary cross jets, attached flow pattern, and large step change in area at the test section inlet (area 
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ratio= 12.02), serve to reduce the swirl velocities to negligible levels in this region. The effect of the 
cross jets on the swirl is the same as that of flow straighteners, the radial jet issuing from the cross jet 
stagnation point serves as a boundary (porous in this case, solid in the case of flow straighteners) across 
which mass flow is inhibited. 

The effect of the attached flow pattern and large step change in area can be explained by angular mo- 
mentum concepts. As the fluid issues from the swirler it possesses a certain amount of swirl, or angular, 
momentum. Since the flow is a wall jet, the fluid moves towards the cell boundaries upon exit from the 
swirler. When the fluid turns in the downstream direction, the large change in area will cause a large 
reduction in the particle velocities due to the law of conservation of angular momentum. These veloci- 
ties may be further reduced due to viscous effects from the cell boundaries, but this would be a small ef- 
fect. This absence of swirl, and the corresponding absence of the decay of swirl, indicates that the recir- 
culation is formed and driven by the step change in area. 

Further downstream of the z=3o0 in. location, at x=4.00 in., the well-defined 4 cell flow pattern is seen 
to deteriorate. Four vortex-type motions are still evident, but the cores of the vortices have translated 
while their strengths have been reduced. The 4.50 in. location shows almost complete breakdown of the 
four cells, with only one very weak region of vortex motion still existing. 

The x=6.00 in. plane illustrates the decay of nearly all motion in the cross-stream plane, the only mo- 
tion being a small amount of flow in the negative z and positive y directions in the lower z range. The 
x=9.00 in. plane shows no detectable motion of any type, this demonstrates the large amount of dissipa- 
tion present in the flowfield as only 6.00 in. upstream existed a very complicated, high speed flow pat- 
tern. Additional line plots of the three mean velocity components are seen in Figures 3.1.1-45 through 
3.1.1-52. 

3. 1.1.3 Turbulence Flowfield 

The turbulence flowfield of the combustor model is best represented by contour plots of the characteris- 
tics in the yz planes for the various downstream locations. The turbulence quantities calculated from 
the individual data samples are the three rms velocities, Unns/ Vrms/ ar *d ^rms> and the correspond- 
ing turbulence kinetic energy, K, as well as the xy component of the Reynolds shear stresses, LTV'. 

The rms velocity and turbulence kinetic energy distributions for the various downstream stations can be 
found in Figures 3.1.1-53 through 3.1.1-99. The data for the x=3.00 in. case consist of only the U rrn s and 
Vrms/ s ‘ nce Wrms was not measurable due to geometric restrictions. 

The profiles for x=1.00 in. can be seen in Figures 3.1.1-53, 3.1.1-54, 3.1.1-55, and 3.1.1-56. The three rms 
velocities seem to exhibit similar behavior with flat contours near the center of the duct. Urms and 
V r ms show a large amount of data scatter for most of the plane, especially the lower z range, while the 
Wrms contour is relatively flat for most of the plane, with scatter only found in the 5.80 in. < z < 6.40 in. 
range. The K distribution, as well, shows a well-defined central region of reduced magnitude extending 
into both corners on the z=9.00 in. boundary. The intensity of the turbulence peaks-up near the walls 
and near the mid y range at the z=9.00 in. boundary. This flat (nearly constant magnitude) central 
region corresponds to the strong region of backflow found at the x=1.00 in. plane, as seen from the mean 
flowfield data, while the peaking near the walls corresponds to the shear layer produced by the strong 
axial flow along the walls. This pattern of reduced levels of turbulence inside the backflow region as 
compared to the shear layer was also seen in the work of Brum and Samuelson (1987) and Vu and 
Gouldin (1982). 

Another characteristic of the K distribution at the x location is the formation of a circular shape, high 
gradient region found in the lower z range. This appears to be the result of the high shear zones from 
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the annular jet flow. This ridge of high K extends around the central zone only for the lower z range and 
does not exist in the upper z range. 

At x=1.50 in. (Figures 3.1.1-57, 3.1.1-58, 3.1.1-59, and 3.1.1-60) the Uj-ms/ Vrms/ an< ^ ^rms contours all 
display relatively flat distributions, resulting in a flat K contour. The circular, high gradient region 
found at x=1.00 in. is barely visible in the lower z range, while a peak is evident near y=2.0 in. on the 
z=9.00 in. boundary. This peak results from both Urms and Vims showing general peaks in that region. 
One characteristic of Vrms illustrated in the contour is the magnitude deficit in the lower right comer. 

Moving to the x=2.00 in. station, some interesting flow characteristics (Figures 3.1.1-61, 3.1.1-62, 

3.1.1- 63, and 3.1.1-64) can be observed. Regarding the line plots, it can be seen that all three rms terms 
(except for Vrm S in the upper z range) and K show very similar behavior across the cell width with 
nearly equal magnitudes (peaking near the mid-y region) and a trailing off of the magnitudes near the 
walls. The velocity gradients are larger in the y direction than in the z direction. This behavior seems 
to be a combination of isotropy, all three rms terms being equal, and homogeneity in the z direction, rms 
terms at a constant value). 

The contours of the turbulence characteristics at x=2.50 in. (Figures 3.1.1-65, 3.1.1-66, 3.1.1-67, and 

3.1.1- 68), show general similarity to those at x=2.00 in., with general peaking in the central region of 
the ceil cross section and larger gradients in the y direction than the z direction. The larger gradients 
result from the close proximity of the measurement plane to the stagnation point; the radial jet fluid is 
concentrated in a thin layer in the y direction while spreading in all direction in the xz plane. 

The behavior of U rm s at x=3.00 in. (Figure 3.1.1-69) does not show much effect from the cross jet injec- 
tion, as most of the contour displays a flat profile with a small peak in the lower right comer. The ef- 
fect of the jets is clearly evident, though, in the Vrms distribution, with a clear and sharp peak occur- 
ring near the center of the cell cross section. Figure 3.1.1-70 presents the contour and line plots of Vnns at 
x=3.00 in. 

Unlike that at the x=3.00 in. location, the Urms contour for x=3.50 in. (Figure 3.1.1-71) shows a large 
central peak with the magnitude dropping off at the cell boundaries. The Vnns distribution given in 
Figure 3.1.1-72 shows magnitudes nearly equal to those at x=3.00 in., with a central peak trailing off at 
the boundaries and walls. The W^s contour, given in Figure 3.1.1-73, however, shows a very asymmet- 
ric behavior, with a general peaking located off center in the negative z direction. Due to this behav- 
ior in W rms , the K distribution also shows an asymmetric peaking in the same general location. Other 
than this abnormality, the distribution shows the turbulence kinetic energy (Figure 3.1.1-74) to be 
highly concentrated in the central region of the duct. 

The behavior of the Urms, V rm s/ W r ms/ and K at the x=4.00 in. station, seen in Figures 3.1.1-75, 

3.1.1- 76, and 3.1.1-77, follows that of the data at x=3.50 in., although with reduced magnitudes due to 
dissipation and decay. The asymmetric Wrms behavior exists in the same general form as at x=3.50 in., 
thus reducing the probability of measurement error. The contour of the turbulence kinetic energy again 
shows the high concentration near the center of the duct, trailing off to nearly constant values at the 
upper and lower walls and showing smooth profiles, peaking at mid y, on the two z boundaries. The 
characteristics at x=4.50 in. (Figures 3.1.1-78, 3.1.1-79, 3.1.1-80, and 3.1.1-81) display the same general 
behavior as the data at both the x=3.50 in. and the x=4.00 in. planes, only with reduced magnitudes. 

The contours of data at x=6.00 in. (Figures 3.1.1-82, 3.1.1-83, 3.1.1-84, and 3.1.1-85) indicate near 
isotropy and homogeneity, although some scatter of the data is evident and the contour does 

show a slight gradient in the negative z direction. The same behavior is seen at x=9.00 in. (Figures 

3.1.1- 86, 3.1.1-87, 3.1.1-88, and 3.1.1-89), with the magnitudes of all three rms terms reduced from those 
at x=6.00 in. Regarding the K contours at both the x=6.00 in. and x=9.00 in. planes, it can be seen that 
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the distributions at x=9.00 in. seem to indicate greater isotropy and homogeneity than at x=6.00 in. 

Since the magnitudes at these stations are so much less than at x=3.50 in., then isotropy and homogene- 
ity can be assumed to be valid assumptions. Although the values of the rms velocities and K at these 
two downstream stations are much lower than the values at the upstream stations, the values are still 
very high when compared to fluid flows such as channel or pipe flow. 

The contour plots of the Reynolds stress component, U r V', can be found in Figures 3.1.1-90 through 
3.1.1-99. Higher levels, either positive or negative, of U'V’ are found in regions of high shear or 
momentum transfer in the xy plane. This component of the Reynolds shear stresses was the only one 
measured. 

From the U’V' contour at x=1.00 in. (Figure 3.1.1-90), the circular ridge of increased U’V' in the left half 
of the plane (looking downstream) can immediately be seen. This ridge was also seen in the K contour 
at x=1.00 in., and corresponds to the high speed flow issuing from the swirlers creating a region of high 
shear and momentum transfer. This ridge appears to indicate that the flow exiting the swirier is not of 
the wall jet variety and appears to possess some type of swirling motion. Near the center of the duct 
one also sees a region of near zero magnitude, also shifted in the positive z direction, signaling a region 
of low or zero shear. This region corresponds to the strong backflow found at this location. 

The x=1.50 in. station exhibits behavior similar to x=1.00 in. only in the sense that a slight ridge is 
visible for the lower z range. In general, the contour is rather flat but asymmetric with a peaking on 
the z=9.00 in. boundary. Also peaking at this same location were Unns and Vnns (and also K), indicat- 
ing the occurrence of some strong turbulence interaction. 

Some interesting behavior appears at x=2.00 in. From Figure 3.1.1-92 regions of large, positive U'V' in 
most of the right half plane and in the lower left comer can be seen, while between these two peaks a 
region of large, negative U'V' can be found. Regarding the xy and yz plane vector plots, the largest 
gradients of U and V can be seen to occur at these locations. At x=2.50 in. the behavior of U'V' is very 
asymmetric, showing a strong positive peak towards the z=9.00 in. boundary. The rest of the U'V’ con- 
tour shows small local peaking but, in general, a flat profile. 

The effect of the cross jets on the turbulent shear can be clearly seen at the x=3.00 in. location. The cen- 
tral portion of the cell cross section can be seen to display large negative values of U'V’, while the 
lower v region shows a trend towards near zero values. The line plot in Figure 3.1.1-94 shows a general 
leveling off near the y boundaries at the y=7.20 in. location, with large slopes occurring near the center 
of the duct. These slopes, and their changes in sign, indicate the strong effect of the cross jet stagnation 
point on the local flowfield. The magnitude of the turbulence stress can be seen to decrease near the cell 
walls and boundaries, where the influence of the cross jets is much reduced. 

The U'V' contour at x=3.50 in. shows two very sharp, positive peaks near the center of the cell cross sec- 
tion, the two sharp positive peaks being separated by a region of much reduced turbulence stress. At the 
x=3.00 in. plane, the positive peak was found in the lower y range, with the much larger negative peak 
occurnng near the center of the duct. For the x=3.50 in. case, the positive peaks are found near the mid- 
dle of the duct with the negative peaks occupying the lower y range. Again the values are seen to trail 
off near the cell walls and boundaries. 

The behavior of the turbulence stress contours at the x=4.00 in. and x=4.50 in. planes is much the same. 
The twin positive peaks occurring at x=3.50 in. are no longer evident as one large peak of nearly the 
same magnitude has developed at x=4.00 in. The location of this peak has translated to the y=1.80 in. 
location, while the negative peak, of much reduced magnitude, remains at the same vertical location. 
The x=4.50 in. location shows a much dissipated peak near y=2.00 in., with the negative peak remain- 
ing at the same location and strength as at x=4.00 in. As with the other cases, the turbulence stress 
drops off near the walls and boundaries. 
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Dissipation and mixing have a large effect on the downstream flowfield, as the x=6.00 in. station 
shows a general spread of the peaks resulting in a flatter U'V' profile. The negative peak is no longer 
evident and, as before, the turbulence stress goes to zero at the walls and boundaries. By x=9.00 in., the 
turbulence stress can be seen to be almost completely dissipated. 

3.1.2 Five-Swirler Water Rig 

The experimental apparatus is shown in Figures 3.1. 2-1, 3. 1.2-2, and 3. 1.2-3. The test section for this 
test is a model of a gas turbine annular combustor. The effects of curvature have been assumed to be neg- 
ligible and the annular flow channel was modeled as a rectangular channel (3 in. x 15 in.) constructed of 
clear plexiglass. At the top of the channel are five annular jets which would provide the fuel/ second- 
ary air mixture in the real combustor. Located 1 channel width, 3 in., downstream of the annular jets are 
a pair of opposing round jets for each annular jet, aligned with the annular jet. These round jets would 
provide the primary air in the actual combustor. The primary jets had dimensions with a diameter 
equaling 0.43 in. and a length/diameter equal to 25. The annular jet had an inner diameter of 1.084 in. 
and an outer diameter of 1.459 in. The annular jets were located 3 in. apart, center to center. The model 
consists of five annular jet/primary jet cells to reduce endwall effects. All measurements were made in 
the center cell. 

The total water flow rate for all jets was determined from the venturi downstream of the test section. 
The pnmary jets' flow rate was determined from the pressure drop in the supply lines to each bank of 
jets. The Ventun and each jet supply line were connected to manometers for flowrate determination. 
Figures 3. 1.2-4 and 3. 1.2-5 present the calibration curves for each flow measurement. The calibration 
data points are shown along with the equations which were used to determine the flow rates from the 
manometer readings. The total flow Venturi output was determined from the difference between the 
manometer readings at the Venturi’s large and small diameters. The primary jet flow pressure drops 
were referenced to the constant head tank level. Readings were taken from the scale mounted next to 
the manometer tubes with divisions that were used as units. 

The concentration field in the model combustor was first measured with the annular jet flow turned off. 
This configuration demonstrated the interactions between two opposing impinging round turbulent jets. 
The total flow rate was 1.7 lb/sec, giving a Reynolds number based on jet diameter of 8000 for each jet. 
The flow control valves were adjusted to position the interface between the jets at approximately the 
center of the channel. Figure 3. 1.2-6 presents the coordinate system used for defining locations within 
the test section. Measurements were made with the laser sheet at z locations of 6.0, 6.5, 7.0, and 7.5 in. 
These locations span the area from the center jets axes to the midway point between adjacent jet pairs. 
Symmetry was assumed and measurements were made on only one side of the jet axes. 

The annular jets were turned on for the second configuration. This allowed visualization of the mixing 
between the primary jets and the annular jet. The total annular jet flow (five jets) was 1.57 Ib/sec and 
the primary jet flow was 1.10 lb/ sec, giving Reynolds numbers based on a hydraulic diameter of 2500 
and 5200 for each annular jet and pnmary jet, respectively. A higher annular jet Reynolds number could 
not be attained without washing the pnmary jet flows downstream and eliminating the interaction be- 
tween the three jets. The ratio of total primary jet momentum to annular jet momentum was 1.26. The 
primary jet flow rate to annular jet flow rate ratio was 0.7. Concentration measurements were made 
with the dye being injected first through one primary jet and then through the annular jet. As with the 
no annular jet flow case, the laser sheet was positioned at z locations of 6.0, 6.5, 7.0, and 7.5 in. for con- 
centration measurements. 

3. 1.2.1 Primary lets Only 

The results from the measurements made with no annular jet flow are presented in Figures 3. 1.2-7 
through 3.1.2-20. At each z location, the data were reduced along lines of constant y, for y values of 0.5, 
1.0, 1.5, and 2.0 in. This provided radial distributions across the jet entering at y = 0.0 in. In addition to 
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these radial profiles, data were reduced along a line of constant x, at x = 3.0 in., corresponding to the jets 
axes and lines parallel to the axes. 

Mean and rms concentration profiles were calculated along each of the lines mentioned above. The 
mean concentration values are presented as a ratio to the maximum mean value on the jet centerline, 

0:1, max- The rms quantities are divided by the local mean value to provide a fluctuation intensity 
quantity, c'/C. 

Figures 3. 1.2-7 through 3.1.2-10 are visualization images of the primary jet concentration field. Figures 

3. 1.2- 7 and 3.1.2-8 are a single frame and a 127 frame average, respectively, with both figures being 
taken through the jet axis. The single frame image is much different from the average image, indicat- 
ing the highly turbulent nature of this flow field. The spread of the jet at the impingement point is 
barely observable in the average image because it switches from top to bottom and in the average is 
greatly reduced. Figures 3.1.2-9 and 3.1.2-10 are 0.5 in. from the jet axis. The spread from the jet im- 
pingement can be clearly seen in the single frame image (Figure 3. 1.2-9). Due to the very turbulent na- 
ture of this interface, the mean image (Figure 3.1.2-10) shows that the concentration is spread over a 
relatively large area. 

Figures 3.1.2-11 through 3.1.2-15 present the mean concentration profiles at the four z locations of the 
laser sheet. In each figure the profiles are presented as line plots and as a 3-D graph. Also included on 
each figure is a diagram indicating where the profiles are located in the test section. The concentration 
profiles in Figure 3.1.2-11 are calculated across the jet axis and are similar to those expected for a single 
jet at y = 0.5 and 1.0 in. At higher values of y, the decrease at the location of the jet axis (x = 3.0 in.) is 
caused by the opposing jet which does not contain dye. With decreasing z, the profiles flatten as the 
dye mixes out toward a constant value. Figure 3.1.2-15 shows the axial mean profile along the jet axis 
and lines parallel to the axis. The jet spreads in the xz plane when it impacts the opposing jet. The 
presence of this spread can be seen in the z = 7.0 in. profile. The mean appears to increase initially from 
the exit value. This effect was also observed by Borleteau (1983) and Long, Chu, and Chang (Tennekes 
and Lumley, 1973). This is probably due to the camera response which is expected to be attenuated at 
the image edge. Long, Chu, and Chang noticed that even after correction of original digital images, the 
response of the camera was still attenuated at the edges. The laser sheet reflected off of the plexiglass 
walls, resulting in a higher laser intensity near the wall. This introduced a nonlinearity near the 
walls which caused the signal to be attenuated. The combination of the camera response and the reflec- 
tion from the plexiglass walls results in the signal attenuation near the walls. 

Figures 3.1.2-16 through 3.1.2-20 present the fluctuation intensity profiles for the case with no annular 
jet flow. Figure 3.1.2-16 is taken through the jet axis and shows a profile similar to that expected for a 
single jet at y = 0.5 and 1.0 in. The interface between the two jets is apparently between y = 1.0 and 1.5 
in. The large peak which would occur at the interface is not observable. The mixing region can be seen 
at y = 1.5 in. for the decreasing values of z. Figure 3.1.2-16 is measured along the jet axis and lines par- 
allel to the axis. The peak, at y = 1.25 in. for the z = 7.5 in. profile, indicates that the jet/ jet interface 
is in fact between y = 1.0 and 1.5 in. The increase for higher values of y is due to the near zero mean 
value which is the denominator of the plotted quantity. 

3. 1.2.2 Annular and Primary lets 

The results for the configuration with annular jet flow and primary jet flow are presented in Figures 

3.1.2- 21 through 3.1.2-51. For the test with dye injection through a primary jet, data reductions were 
performed along radial and axial lines as described for the no annular jet flow case. The test with dye 
in the annular jet flow was analyzed by reducing data along lines of constant x, for x values of 1.0, 2.0, 
3.0, and 4.0 in. This provided radial distributions across the annular jet. An axial profile was also cal- 
culated along the annular jet's axis (y = 1.5 in.) and lines parallel to the axis. Mean values were nor- 
malized by the maximum annular jet centerline mean value, Cel, max- 
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Figures 3.1.2-21 through 3.1.2-35 present the results from the test with the dye being injected into a pri- 
mary jet. Figures 3.1.2-21 through 3.12-24 are visualization images of the primary jet concentration 
field. As with the previous case, the average images are greatly different from the single images indi- 
cating a highly unsteady flowfield. The effect of the annular jet is dearly evident in these images. 

The primary jet flow is driven downstream by the annular jet. The figures at z = 7.0 in. show much less 
primary jet flow at this location than in the case without annular jet flow. This indicates that the 
primary jet is more thoroughly mixed at z = 7.0 in., due to the annular jet flow. 

Figures 3.1.2-25 through 3.1.2-29 present the mean concentration profiles. Figure 3.1.2-25 shows the pro- 
files through the jet axis. These are similar to the case without annular jet flow, except that above the 
primary jets the mean value is nearly zero due to the annular jet flow. The mean values drop off with 
decreasing z as before. This case results in a faster decrease due to the mixing with the annular jet. 
Figure 3.1.2-29 shows profiles along the jet axis and lines parallel to the axis. The value at the edge is 
attenuated as discussed previously. 

Figures 3.1.2-30 through 3.1.2-35 present the fluctuation intensity profiles. Figure 3.1.2-30 shows pro- 
files similar to the expected curves for a single jet at y = 0.5 and 1.0 in. The peak on the lower x side of 
the jet is higher due to the increased mixing with the annular jet. This can be still be seen at z = 7.0 in. 
(Figure 3.1.2-31). At z = 6.5 in., the fluctuation intensity profiles are nearly zero, indicating that the 
mixing is nearly zero. Figures 3.1.2-34 and 3.1.2-35 are both fluctuation intensity profiles along the 
plane of the jet axis. Figure 3.1.2-34 only goes to y= 1.7 in. because the quantity c'/C was very high and 
erratic beyond that point, due to the nearly zero mean value. Figure 3.1.2-35 shows the rms normalized 
by Cd, max instead of C to give a better indication of the shape of the rms profile. 

Figures 3.1.2-36 and 3.1.2-37 compare the primary jet concentration results from the test cases with, and 
without, annular jet flow. Figure 3.1.2-36 compares mean and fluctuation intensity profiles across the 
primary jet axis (z = 7.5 in.). The effect of the annular jet flow is seen in the mean as a decrease on the 
low x value side of the jet. The fluctuation intensity peak on this side of the jet is higher, with the an- 
nular jet flow, due to the increased mixing. Figure 3.1.2-37 presents mean and fluctuation intensity pro- 
files along the primary jet axis and lines parallel to the axis. The mean profiles at z = 7.5 in. show that 
the annular jet flow does not have a very large effect at this location. That is because the annular jet 
momentum, although greater, is spread over a larger area than that of the primary jet. The fluctuation 
intensity peak at the intersection of the three jets is higher due to a more turbulent flow with the 
annular jet flow. At the lower z values, the profiles for the case with annular jet flow lower mean val- 
ues and higher fluctuation intensity values. Both of these trends indicate more mixing with the annu- 
lar jet flow, as expected. 

The results for the case with dye in the annular jet flow are presented in Figures 3.1.2-38 through 3.1.2- 
51. Figures 3.1.2-38 through 3.1.2-41 are concentration field visualization images of the annular jet. 

The unsteady nature of this field can be deduced from the large differences between the single frame 
image and the average image. This flow does not appear as turbulent as the primary jet flow and in- 
deed the Reynolds number is significantly lower for the annular jet flow. The images appear asymmet- 
ric with more dye on the right side of the flow. This is apparently due to an uneven flow distribution at 
the annular jet exit. This could be owing to the annular jet being nearly completely closed. If the jet con- 
trol valve had been opened wider, the higher annular jet flow would have washed the primary jet 
flows downstream and eliminated any interaction between the jets in the field of observation. 

The mean concentration profiles are presented in Figures 3.1.2-42 through 3.1.2-146. Figure 3.1.2-42 
shows the mean profiles through the annular jet axis. The double peaked curve at x = 1.0 in. is slightly 
asymmetric, indicating the uneven flow distribution discussed above. The profile flattens and spreads 
downstream until, downstream of the primary jets at x = 4.0 in., the profile is essentially flat. At z = 7.0 
in.,= there is only a single peak because the laser sheet intersects the annulus of the annular jet in just 
one area. This peak is off center due to the uneven flow distribution from the annular jet. Again, the 
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profile spreads and flattens downstream. At z = 6.5 in., the profile increases at x = 3.0 in. due to the 
primary jets carrying the annular jet flow across the channel to this location (Figure 3.1.2-44). The in- 
crease is not as significant at z = 6.0 in. due to the increased distance from the primary jet pair and the 
effect of the adjacent primary jet pair; this location is midway between two jet pairs. Figure 3.1 .2-46 
shows mean profiles along the annular jet axis and lines parallel to the axis. The effect of the primary 
jet pair can be seen in the profiles at z = 7.0 and 6.5 in. as an increased concentration level extending 
across the channel at x = 3.0 in. 

Fluctuation intensity profiles are presented in Figures 3.1.2-47 through 3.12-51. Figure 3.1.2-47 presents 
measurements made with the laser sheet aligned with the annular jet axis. The edges of the jet can be 
seen as peaks. The overall level increases at x = 3.0 in., due to the mixing with the primary jets at this 
location. The level drops again at x = 4.0 in., downstream of the primary jets. The fluctuation intensity 
levels decrease with decreasing z as the concentration reaches a steady value and the mixing drops off. 
At z = 6.0 in., the curves are low in the middle and increase towards the sides due to the low mean value 
in these areas. Figure 3.1.2-51 presents fluctuation intensity profiles along the annular jet axis and lines 
parallel to the axis. The profiles at z = 7.5 in., the annular jet axis, have the expected shape, nearly 
zero until the primary jet location at x = 3.0 in. The peak at this location indicates the highly turbulent 
mixing occurred where all these jets interacted. 
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Figure 3.1 .1-1. Three-view drawing of annular combustion chamber model. 
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Figure 3.1 .1-2. Photograph of annular combustor chamber model. 



Figure 3. 1.1-3. Arrangement of the swirler cells. 
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Figure 3. 1.1-4. Actual annular combustor. 
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Figure 3. 1.1-5. Allison Gas Turbine 570- K annular combustor. 



Figure 3. 1.1-6. Annular jet and throttle valve. 


49 




n 


u 

TE92-2259 


Figure 3.1. 1-7. Vane swirier from model combustor. 
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Figure 3.1. 1-8. Original swirier orientation and flow pattern. 
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Figure 3. 1.1-9. Modified swirler orientation and flow pattern. 



Figure 3.1.1-10. Primary cross jet and throttle valve. 
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Figure 3.1.1-11. Plenum chamber and centrifugal fan. 
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Figure 3.1.1-12. Contour and line plots of V for primary cross jet at y=0.06 in. 
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Figure 3.1.1-13. Contour and line plots of V^ns for primary cross jet at y=0.06 in. 
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Figure 3.1.1-14. Comparison of V data at y=0.06 in. to classical pipe flow data of Laufer. 
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Figure 3.1.1-16. Contour and line plots of V for primary cross jet at y=0.25 in. 
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Figure 3.1.1-17. Contour and line plots of Vnns for primary cross jet at y=0.25 in. 
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Figure 3.1.1-18. Contour and line plots of V for primary cross jet at y=0.50 in. 
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Figure 3.1.1-19. Contour and line plots of Vfms for primary cross jet at y=0.50 in. 
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Figure 3.1.1-20. Contour and line plots of V for primary cross jet at y=0.75 in. 
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Figure 3.1.1-21. Contour and line plots of Vrms for primary cross jet at y=0.75 in. 
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Figure 3.1.1-22. Contour and line plots of V for primary cross jet at y=l. 00 in. 
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Figure 3.1.1-23. Contour and line plots of Vrms for primary cross jet at y=1.00 in. 
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Figure 3.1.1-24. Contour and line plots of V for primary cross jet at y=1.25 in. 
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Figure 3.1.1-25. Contour and line plots of Vnns for primary cross jet at y=1.25 in. 
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Figure 3.1.1-26. Contour and line plots of V for primary cross jet at y=1.50 in. 
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Figure 3.1.1-27. Contour and line plots of Vnm for. primary cross jet at y=1.50 in. 
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Figure 3.1 .1-28. xy plane measurement grid for the primary cross jets. 
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Figure 3.1.1-29. xz plane measurement grids for the primary cross jels for: (a) y=0.06 in. and y=0.25 in. 

(b) y =0.50, 0.75, 1.0, 1.2. 
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Figure 3.1.1-30. yz plane measurement grid. 
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Figure 3.1.1-31. xy plane measurement grid. 
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Figure 3.1.1-32. xy plane mean velocity vector plots: (a) z=5.80 in. (b) z=6.00 in. (c) z=6.20 in. 
(z = 6.0 in. is the midplane between two swirlers). 
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Figure 3.1.1-33. 


xy plane mean velocity vector plots: (a) z=6.40 in. (b) z=6.60 in. (c) z=6.80 in. 
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Figure 3.1.1-34. xy plane mean velocity vector plots: (a) z=7.00 in. (b) z=7.20 in. (c) z-7AQ in. 
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Figure 3.1.1-35. xy plane mean velocity vector plots: (a) z=7.60 in. (b) z=7.80 in. (c) z=8.00 in. 
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Figure 3.1.1-36. xy plane mean velocity vector plots: (a) z=8.20 in. (b) z=8.40 in. (c) z=8.60 in. 
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Figure 3.1.1-37. xy plane mean velocity vector plots: (a) z=8.80 in. (b) z-9.00 in. 
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Figure 3.1.1-38. Recirculation zone cross-sections: (a) x=1.00 in. (b) x=1.50 in. (c) x=2.00 in. (d) x=2.50 

in. 
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Figure 3.1.1-39. y v. plane mean velocity vector plots: (a) x=l.(X) in. (b) x=1.50 in. 




/ 

* 

\ 

\ 


•. V V - 

» V \ 

.. V x 

^ \ r . 

. T / 

^ / X . 


/ 

/ l 
• ‘ V 

' \ 

' ' \ 
1 1 V 

t 


\ 

\ 

\ 

l 

s 


l 

l 

l 

•» 

\ 


\ ^ ' 
M I*'' '' 


•- V 




*7 


> 

i V 


\ 

\ 

t 

l 


i l l 

l v - 


r j 


o 

in 


CM 

05 

CM 

CM 

CM 

05 




' ' S S 


/ 

/ 


/ / / 


/ / / ' 
/ / / <- 

‘ ' L ' 

» \ \ ’ 

\ - 
• 5 » - 


• ‘ \ * 

’ll' 

< / l 5 

'‘I s 
' ‘ l ‘ 


„ / / 

j j f 


v < i 


/ J- » 



83 


Figure 3.1.1-40. y z plane mean velocity vector plots: (a) x=2.00 in. (b) x=2.50 in. 
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Figure 3.1.1-41. YZ plane mean velocity vector plots: (a) x=3.5() in. (b) x=4.00 in. 






•*r 

<y> 

CNJ 

CM 

CM 

Ol 

UJ 




I i i 

• * + s I i i t i 

'// / / // y 

‘ ' Jill'' ' 

* / l l t ' ' ' 

■/////''' 
'/////'* ' 

i. t ‘ • • 


I 



o 

m 


85 


Figure 3.1.1-42. yz plane mean velocity vector plots: (a) x=4.50 in. (b) x=6.00 in. 
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Figure 3.1.1-43. yz plane mean velocity vector plot at x=9.00 in. 
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Figure 3.1.1-44. Radial jet flow pattern: (a) xy plane view (b) yz plane view. 
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Figure 3.1.1-45. Line plo! of U at X=100, 1.50, 2.00, and 2.50 in. 
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Figure 3.1.1-46. Line plots of U at x=3.00, 3.50, 4.00, and 4.50 in. 
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Figure 3.1.1-47. Line plots of U at x=6.00 in. and 9.00 in., and of V at x=1.00 in. and 1.50 in. 
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Figure 3.1.1-48. Line plots of V at x=2.00, 2.50, 3.(X), and 3.50 in. 
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Figure 3.1.1-49. Line plots of V at x=4.00, 4.50, 6.00, and 9.00 in. 
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Figure 3.1.1-50. Line plots of W al x=l.Q0, 1.50, 2.00, and 2.50 in. 
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Figure 3.1.1-51. Line plots of W at x=3.50, 4.00, 4.50, and 6.00 in. 
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Figure 3.1.1-52. Line plots of W at x=9.00 in. 
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Figure 3.1.1-53. Contour and line plots of Unns at x=1.00 in. 
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Figure 3.1.1-54. Contour and line plots of Vnns at x=1.00 in. 


97 





Wrms (ftVsec.) 



-& 



Y (in.) 


TE92-2307 


Figure 3.1.1-55. Contour and line plots of at x=1.00 in. 
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Figure 3.1.1-56. Contour and line plots of K at x=1.00 in. 
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Figure 3.1.1-57. Contour and line plots of Unns at x=1.50 in. 
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Figure 3.1.1-58. Contour and line plots of Vnj,; at x=1.50 in. 
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Figure 3.1.1-59. Contour and line plots of W-rms at x=150 in. 
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Figure 3.1.1-60. Contour and line plots of K at x=1.50 in. 
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Figure 3.1.1-61. Contour and line plots of Urms at x=2.00 in. 
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Figure 3.1.1-62. Contour and line plots of Vnns at x=2.00 in. 
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Figure 3.1.1-63. Contour and line plots of Wnns at x=2.00 in. 
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Figure 3.1.1-64. Contour and line plots of K at x=2.00 in. 
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Figure 3.1.1-65. Contour and line plots of Unns at x=2.50 in. 
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Figure 3.1.1-66. Contour and line plots of Vnrts at x=2.50 in. 
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Figure 3.1.1-67. Contour and line plots of Wnr\s at x=2.50 in. 
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Figure 3.1.1-68. Contour and line plots of K at x=2.50 in. 
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Figure 3.1.1-69. Contour and line plots of Urms at x=3.00 in. 


112 




Vrms (fL/sec.) 




Y (in.) 

TE92-2322 


Figure 3.1.1-70. Contour and line plots of Vrms at x=3.00 in. 
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Figure 3.1.1-71. Contour and line plots of Unrvs at x=3.50 in. 
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Figure 3.1.1-72. Contour and line plots of Vj^g at x=3.50 in. 
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Figure 3.1.1-73. Contour and line plots of Wnns at x=350 in. 
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Figure 3.1.1-74. Contour and line plots of K at x=3.50 in. 
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Figure 3.1.1-75. Contour and line plots of Urms at x=4.00 in. 
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Figure 3.1.1-76. Contour and line plots of Vrmg at x=4.00 in. 
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Figure 3.1.1-77. Contour and line plots of Wnus at x=4.00 in. 
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Figure 3.1.1-78. Contour and line plots of K at x=4.00 in. 
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Figure 3.1.1-79. Contour and line plots of Urms at x=4.50 in. 
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Figure 3.1.1-80. Contour and line plots of Vnns at x=4.50 in. 
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Figure 3.1.1-81. Contour and line plots of K at x=4.50 in. 
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Figure 3.1.1-82. Contour and line plots of Urms at x=6.00 in. 
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Figure 3.1.1-83. Contour and line plots of Vrms at x=6.00 in. 
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Figure 3.1.1-84. Contour and line plots of Wnus at x=6.00 in. 
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Figure 3.1.1-85. Contour and line plots of K at x=6.00 in. 
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Figure 3.1.1-86. Contour and line plots of Urms at x=9.00 in. 
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Figure 3.1.1-87. Contour and line plots of Vnns at x=9.00 in. 
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Figure 3.1.1-88. Contour and line plots of Wrms at x=9.00 in. 
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Figure 3.1.1-89. Contour and line plots of K at x=9.00 in. 
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Figure 3.1.1-90. Contour and line plots of U'V' at x=1.00 in. 
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Figure 3.1.1-91. Contour and line plots of U’V' at x=150 in. 
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Figure 3.1.1-92. Contour and line plots of U’V at x=2.00 in. 
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Figure 3.1.1-93. Contour and line plots of U'V’ at x=2.50 in. 
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Figure 3.1.1-94. Contour and line plots of U'V at x=3.00 in. 
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Figure 3.1.1-95. Contour and line plots of U’V' at x=350 in. 
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Figure 3.1.1-96. Contour and line plots of U'V’ at x=4.00 in. 
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Figure 3.1.1-97. Contour and line plots of U'V’ at x=4-50 in. 
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Figure 3.1.1-98. Contour and line plots of U'V at x=6.00 in. 
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Figure 3.1.1-99. Contour and line plots of U'V at x=9.00 in. 
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Figure 3. 1.2-1. Model annular combustor test rig schematic. 
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Figure 3.1 .2-2. Detail of model annular combustor test section. 
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Figure 3.1 .2.-3. Section of model annular combustor test section. 
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Figure 3.1. 2-4. Flow calibration curves for primary jet flow. 
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Figure 3.1.2-5. Venturi calibration curve for total jet flow. 
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Figure 3.1 .2-8. One hundred twenty-seven frame average primary jet concentration field visualization 

image without annular jet flow (z = 7.5 in.). 


Figure 3. 1.2-7. Single frame primary jet concentration field visualization image without annular jet 

flow (z = 7.5 in.). 
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Figure 3.1. 2-9. Single frame primary jet concentration field visualization image without annular jet 

flow (z = 7.0 in.). 



TE92-2362 

Figure 3.1.2-10. One hundred twenty-seven frame average primary jet concentration field visualization 

image without annular jet flow (z = 7.0 in.). 
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Figure 3.1.2-11. Mean primary jet concentration without annular jet flow (z = 7.5 in.). 
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Figure 3.1.2-12. Mean primary jet concentration without annular jet flow (z = 7.0 in.). 
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Figure 3.1.2-13. Mean primary jet concentration without annular jet flow (z = 6.5 in.). 
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Figure 3.1.2-14. Mean primary jet concentration without annular jet flow (z = 6.0 in.). 
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Figure 3.1.2-15. Mean primary jet concentration without annular jet flow (x = 3.0 in.). 
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Figure 3.1.2-16. Root mean square primary jet concentration without annular jet flow (z = 7.5 in.). 
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Figure 3.1.2-17. Root mean square primary jet concentration without annular jet flow (z = 7.0 in.). 
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Figure 3.1.2-18. Root mean square primary jet concentration without annular jet flow (z = 6.5 in.). 
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Figure 3.1.2-19. Root mean square primary jet concentration without annular jet flow ( z = 6.0 in.). 
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Figure 3.1.2-20. Root mean square primary jet concentration without annular jet flow (x = 3.0 in.). 
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Figure 3.1.2-21. Single frame primary jet concentration field visualization image with annular jet flow 

(Z = 7.5 in.). 
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Figure 3.1.2-22. One hundred twenty-seven frame average primary jet concentration field visualization 

image with annular jet flow (z = 7.5 in.). 
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Figure 3.1.2-23. Single frame primary jet concentration field visualization image with annular jet flow 

(z = 7.0 in.). 
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Figure 3.1.2-24. One hundred twenty-seven frame average primary jet concentration field visualization 

image with annular jet flow (z = 7.0 in.). 
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Figure 3.1.2-25. Mean primary jet concentration with annular jet flow (.z = 7.5 in.). 
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Figure 3.1.2-26. Mean primary jet concentration with annular jet flow (z = 7.0 in.). 
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Figure 3.1.2-27. Mean primary jet concentration with annular jet flow (z = 6-5 in.). 
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Figure 3.1.2-28. Mean primary jet concentration with annular jet flow (z = 6.0 in.). 
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Figure 3.1.2-29. Mean primary jet concentration with annular jet flow (x = 3.0 in.). 
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Figure 3.1.2-30. Root mean square primary jet concentration with annular jet flow (z = 7.5 in.). 
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Figure 3.1.2-31. Root mean square primary jet concentration with annular jet flow (z = 7.0 in.). 
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Figure 3.1.2-32. Root mean square primary jet concentration with annular jet flow (z = 6.5 in.). 
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Figure 3.1.2-33. Root mean square primary jet concentration with annular jet flow (z = 6.0 in.). 


175 




LOCATION (INCHES) 


Laser sheet location 

Concentration profile 
locations 



Figure 3.1.2-34. Root mean square primary jet concentration with annular jet flow (x = 3.0 in.). 
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Figure 3.12-35. Root mean square primary jet concentration normalized by C cl. max with annular jet 

flow (x = 3.0 in.). 
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Figure 3.1.2-36. Comparison of primary jet concentration with and without annular jet flow (z = 7.5 in.). 
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Figure 3.1.2-37. Comparison of primary jet concentration with and without annular jet flow ( x = 3.0 in.). 
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TE92-2390 

Figure 3.1.2-38. Single frame annular jet concentration field visualization image (z = 7.5 in.). 
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Figure 3.1.2-39. Frame average annular jet concentration field visualization image (z = 7.5 in.). 
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Figure 3.1.2-40. Single frame annular jet concentration field visualization image (z = 7.0 in.). 
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Figure 3.1.2-41. Frame average annular jet concentration field visualization image (z = 7.0 in.). 
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Figure 3.1.2-42. Mean annular jet concentration (z = 7.5 in.). 
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Figure 3.12-43. Mean annular jet concentration (z = 7.0 in.). 
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Figure 3.1.2-44. Mean annular jet concentration (z = 6.5 in.). 
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Figure 3.12^45. Mean annular jet concentration (z = 6.0 in.). 
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Figure 3.1.2-46. Mean annular jet concentration (y = 1.5 in.). 
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Figure 3.1.2-47. Root mean square annular jet concentration (z = 7.5 in.). 
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Figure 3.1.2-48. Root mean square annular jet concentration (z = 7.0 in.). 
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Figure 3.1.2-49. Root mean square annular jet concentration (z = 6.5 in.). 
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Figure 3.12-50. Root mean square annular jet concentration (z = 6.0 in.). 
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Figure 3.12-51. Root mean square annular jet concentration V y = 1-5 in.). 
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3.2 FLOW VISUALIZATION STUDY 


An extensive flow visualization study was conducted on both water and airflow rigs to delineate the 
turbulent flow structure of the flowfield established by the dome swirler, primary jets, and their inter- 
action. This effort included a systematic investigation of the dynamic flow structure on the five- 
swirler water rig, swirling jet established by a rotating pipe, single swirler, dome annular jets without 
swirl, primary jets, and interaction of primary jets with dome annular jets. 

The turbulent flow structure was quite complex and exhibited considerable time unsteadiness generally 
stipulated for gas turbine combustor flowfield. In addition to the video tape records, a number of still 
flow visualization photographs were taken. The salient features of the flowfield are discussed in ap- 
propriate subsections of Sections III and IV. 


195 


3.3 RIG REFINEMENT ACTIVITIES 


Based on the results from the air and water rigs, two additional experiments on a swirling jet and a sin- 
gle swirler apparatus were performed. 

3.3.1 Swirling let 

Of all the characteristics of swirling flows the recirculation zone is of the most interest. The recircula- 
tion zone is critical in the stability of the combustion process and, hence, must be fully understood in or- 
der to allow the improved design of combustion chambers. Although the size, shape, and strength of 
the recirculation zone depend on the method of swirl application or generation, the general character- 
istics remain the same. 

This segment of the research involved the study of seed particle concentration fields of a swirling, tur- 
bulent jet. Using a light scattering concentration measurement technique known as marker nephelome- 
try, color enhanced video images of these fields were produced. These images allowed the analysis of 
the formation of the recirculation zone and associated characteristics of the swirling jet. Discussions of 
the experimental apparatus along with the presentation and discussion of the results are included in 
this subsection. 

3 . 3. 1.1 Experimental Apparatus 

The experimental jet apparatus for the swirling jet concentration studies is shown in Figure 3.3. 1-1. Am- 
bient laboratory air was drawn into the mixing section by an electric fan, with the injection of the seed 
particles occurring immediately downstream. The seed particles were created by atomization of a com- 
bination of ethyl alcohol and fog juice (Deod. Apco 467) with a TSI Atomizer and entered the mixing 
section of the jet apparatus with an approximate particle diameter of 1 pm. This type of seed particle 
was favored over other possible seed particles due to very good light scattering properties as well as 
the ability of the particles to follow the actual fluid motion. The mixing section contained several 
banks of screen and straws in order to reduce the level of turbulence, the straws and screens reduced the 
turbulence scales and thus increased dissipation. A contraction section (contraction ratio= 36.0) was 
mounted to the exit of the mixing section to increase the axial velocity and further decrease the turbu- 
lence. 

Upon exit from the contraction section the flow entersed the rotating section, a 3-D view of which can be 
seen in Figure 3.3.1-2. The rotating section consisted of a straight tube with an inside diameter of 2.00 
in., the entire length of which was filled with straws. These straws, 3.00 mm in diameter, drove the 
solid body rotation of the flow. The rotation in this section was produced by a hydraulic motor via gear 
and chain drive, allowing infinite adjustability of the rotational speed. The axial average could also 
be controlled, allowing nearly any combination of axial speed, rotational speed, and swirl ratio. The 
mass-averaged speed through the jet chosen for this experiment was 8.40 ft/ sec, giving a Reynolds num- 
ber based on jet diameter of 8900. 

3.3. 1.2 Results And Discussion 

Two-dimensional concentration fields/flow visualization images of the initial region (x/D < 7) of a 
swirling turbulent jet are given for seven different swirl ratios: 0.00, 0.33, 0.66, 1.00, 1.33, 1.66, and 2.00. 
For each of these swirl ratios color enhanced video images corresponding to single frames and 127 aver- 
aged frames are shown. These images can be seen in Figures 3.3.1-2 through 3.3.1-8. 

The imaged region of the flowfield extended roughly 7 jet diameters in both the downstream direction 
and direction perpendicular to the jet axis. With a jet diameter of 2.0 in., the imaged area was 252 in. 
For this set of images the laser sheet thickness was 0.04 in., resulting in a measurement volume size of 
0.03 in. x 0.03 in. x 0.04 in., small enough to resolve all but the smallest scales found in the flowfield, 
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such as those corresponding to the Kolmogoroff universal length scale. Using an estimate of the turbu- 
lent dissipation from Tennekes and Lumley (1973), the Kolmogoroff universal length and time scales 


il = 


V Yt 

V e / 


( 2 ) 


T = 



(3) 


are found to be on the order of 0.01 in. and 0.005 sec. Since the camera operated at a 60 Hz rate, resulting 
in a time resolution of 0.0167 sec., the single frame images are not instantaneous but averaged over a 
small time interval. 


When these images were acquired great care was taken to ensure that the laser intensity was the same 
for all the cases, thereby allowing indirect comparison between images. Before a new case was to be 
imaged, the test facility was flushed with fresh ambient air to ensure a near-zero background marker 
particle concentration level. This was required since the during the operation of the swirling jet appa- 
ratus, the room became rather quickly saturated with the marker particles. The saturation occurred 
very slowly, though, compared with the time required for a complete sequence of images at one swirl 
ratio, allowing characteristics such as the time-averaged spread angle to remain constant in time. 

The variations in the flowfield as the swirl ratio is increased as well as the differences between the 
single frames and multiple frame averages can immediately be seen. In general, several phenomena oc- 
cur with an increase in swirl: 


• an increase in the spread angle of the jet 

• the formation of a recirculation zone when a critical swirl ratio is reached 

• the growth of this recirculation region with an increase in swirl 

The zero swirl case can be seen in Figure 33.1-2. Jet flows with zero swirl are characterized by a very 
small jet spread angle which, for the 127 frame average image, is measured to be 24 deg. From Figure 
3.3. 1-2 (b) some apparently wavy behavior of the jet boundary can be seen. This behavior can be at- 
tributed to the use of straws in the rotating section. Some of the straws were not exactly aligned paral- 
lel to the jet axis and interacted with neighboring straws to develop undesirable secondary flows and jet 
growth. This effect does cause some small errors in the spread angle measurement for the swirl number, 
N=0.00 and N=0.33, cases (swirl number is the ratio of angular momentum to axial momentum), but not 
for the larger levels of swirl. 

The spread angles given for these images represent the spread of the visible radial extent of marker 
particles with downstream distance. Most jet studies, usually involving the measurement of velocity 
profiles, define the spread angle as the spread of the velocity half radius with downstream distance. 

A comparison of the single and multiple frame averaged images reveals the dissimilar nature of the 
single frame and averaged profiles. The single frame image can be seen to show very large-scale fluctu- 
ations for x/D > 2, with a correspondingly large amount of intermittency near the jet boundary, while 
the averaged image shows very smooth contours. Dahm and Dimotakis (1985) discovered the same be- 
havior in their jet studies. 

in Figure 3.3. 1-2 a discrepancy can be seen in the background intensity level, above and below the jet, 
and an apparent distortion of the comers of the images is also visible. The distortion effect is due to 
nonlinear camera response, while the variation in the background intensity level is due to some type of 
reflection causing undesired illumination. These two effects are present in all the processed images but 
do not hinder the study of the flow structure. 
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Figure 3.3.1-3 shows single and averaged images corresponding to the first case of applied swirl, 

N=0.33. The measured spread angle, 22 deg, is smaller than that reported for the zero swirl case but, as 
stated before, the results of these two cases are in slight error. As with the zero swirl case, very large- 
scale fluctuations can be seen from the single frame image, resulting in profiles that do not resemble the 
mean. 

The images for the N=0.66 case can be found in Figure 3.3.1-4. Immediately visible in the averaged im- 
age is the reduced downstream extension of the various intensity levels compared to the N=0-33 and 
N=0.00 cases. This is due to the reduction in the rate of the downstream convection which, in turn, is 
due to the development of a small, axial pressure gradient with increase in swirl. This pressure gradi- 
ent also existed at the N=033 case, but was so small as to be almost negligible. Unlike the N=0.00 and 
N=0.33 cases, which show a nearly constant rate of jet spread with downstream distance, the spread for 
N=0.66 displays two distinct spread angles or 30 deg for x/D < 2.5 and 12 deg for x/D > 2.5. 

The first sign of recirculating flow can be seen in the single frame image for N=1.00 (Figure 3.3.1-5), 
where large-scale fluctuations can be seen to occur much closer to the jet exit than for the previous cases. 
The regions of higher concentration appear to diverge at a downstream distance of x/D=1.5, indicating 
a backflow of fluid with lower marker particle concentration. This fluid must recirculate from stations 
located farther downstream. The averaged image shows no direct signs of recirculation, only a further 
reduction in the downstream extension of the various concentration levels. The averaged image also 
shows the change in the spread angle at the downstream location x/D=2, similar to that at N=0.66. 

The spread angle measured in the initial region of the flowfield is 45 deg. 

The next level of swirl, N=1.33, shows signs of the recirculation zone in both the single frame and the 
multiple frame average images (Figure 3.3.1-6). The fluid exiting the jet, with very high levels of 
marker particle concentration, can be seen to move radially around a region of marker particle deficit, 
which is due to backflow from downstream stations. The formation of the forward stagnation point can 
be clearly seen, but is not distinctly visible for the rear (downstream) stagnation point. As with the 
previous cases, large-scale fluctuations are prevalent in most of the flowfield. An increase in the size of 
the region of intermittency also occurs as small volumes of ambient fluid are transported across most of 
the radial extent of the flowfield. Similar to the previous cases, the spread angle can be seen to have 
an increased value, 80 deg, in the initial portion of the flowfield (x/D < 1), with the subsequent motion 
of the highly concentrated fluid around the recirculation zone. 

The single frame and multiple frame averaged images for the N=1.66 case can be seen in Figure 3.3.1-7. 
The recirculation zone formed at this level of swirl is larger than for the N=1.33 case, with a corre- 
sponding increase in the spread angle (to 140 deg) in the initial mixing region. The single frame image 
shows fluid of very low concentration fluid recirculating. Small regions of unmixed, ambient fluid can 
be seen to exist in the recirculation zone, but at a very low rate of occurrence. As with the N=1.33 case, 
the location of the forward stagnation point can be clearly seen while the rear (downstream) stagnation 
point cannot be discerned. 

Figure 3.3. 1-8 illustrates the concentration field for the N=2.00 case. From this figure the fluid motion 
appears to be completely upstream, as the highly concentrated fluid can be seen to issue from the jet exit 
and flow around the lip of the jet. In actuality, the recirculation zone is so large that the imaged area 
is only a small portion of the total recirculation zone size. This same type of behavior was observed by 
Gore and Ranz (1964), who applied enough swirl to a jet flow to cause it to diverge at right angles to the 
jet axis at the jet exit. The figure shows that most of the fluid existing in the recirculation zone is ambi- 
ent, or of very low concentration levels. This indicates the large amount of diffusion that has occurred 
during the fluid motion around the recirculation zone, another sign of the large scale of the recirculation 
zone. 
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Finally, Figure 3.3.1-9 shows the variation in the spread angle as a function of the swirl ratio. The 
very smooth resulting curve, which increases past 180 deg for the higher swirl ratios, is visible. This 
set of data can be curve fit to be represented by the relation: 

e = 12.49e 141N (4) 

3.3.2 Single Swirler Configuration 

In the combustor it is desirable that the annular swirler produce a strong recirculation zone. This area 
of reverse flow holds the hot gases in the center of the chamber away from the liner. It also produces 
large amounts of fuel/air mixing near the swirler exit and shortens the needed combustor length. 

Each swirler consists of an actual swirler from the production of the Allison 570-K turboshaft gas tur- 
bine engine (Figure 33.2-1). These are 60 deg angle flat vane swirlers with 12 vanes. They have an 
outer flow diameter of 1.459 in. and inner flow diameter of 1.084 in. For nonswirling flow (S = 0), the 
normal swirler is replaced with a dummy swirler having no vanes. 

A measure of the strength of the swirling flow produced is the swirl number. The most commonly ac- 
cepted definition: 


Swirl Number = 


C. 

G x R 2 


(5) 


where: 


Gq angular momentum flux 
Gx axial momentum flux 
R2 outer radius of swirl generator 

For vane swirlers, the swirl number can be approximated by (Beer and Chigier, 1983; Mathur and Mc- 
callum, 1967): 


S 



v R 2 ; 
< R 2 , 


tan0 


where: 


Rl = inner flow radius 
R2 = outer flow radius 
9 = swirler/ vane angle 


( 6 ) 


From the geometry of the swirlers, and using Equation 6, the swirl number is calculated to be 1 .52. Swirl 
generators producing swirl numbers greater than 0.6 produce recirculation regions along the swirling axis 
near the exit (Beer and Chigier, 1983; Gupta et al, 1984). It is expected that the swirlers used will pro- 
duce a very strong recirculation region. The annular swirler inlet tubes are made of stainless steel and 
have diameters equivalent to the flow diameters of the swirlers. Each inlet tube is 25 effective diame- 
ters, 9.547 in. long. Where the effective diameter is defined by: 


D e = outer diameter - inner diameter 


(7) 
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Initially the swirlers were mounted at the end of the inlet annulus flush with the headplate. How- 
ever, in doing some preliminary smoke flow visualization an unexpected flow phenomenon occurred. 

The flow would leave the swirl er and make a sharp 90 deg bend and hug the headplate all the way out 
to the side walls. At the side walls, the flow would then proceed downstream. This phenomenon has 
also been observed by Gore and Ranz (1964) when high degrees of swirl are generated flush with the 
orifice exit. 

It was not believed that this type of flow profile accurately modeled the intended combustor flow. 
Therefore, the swirlers were recessed 0.25 in. back from the headplate. The annulus with the recessed 
swirler is shown in Figure 3.3.2-2. 

The primary jets provide cool additional air for mixing and combustion with the swirling air. Each 
primary jet is made of stainless steel and is 25 diameters long, with the diameters varying for the dif- 
ferent configurations. The necessary jet diameters were determined using the required mass flow ratios 
and number of jets. Referring to Figure 3.3.2-3, the mass flow ratio is given by: 

Mj N iVj(^k (g 

M s N s v s (A t N b )/N b 

Imposing the condition that Vj is equal to V s and solving 3.8 for dj yields: 




4MjN s A, 

reM s Nj 


(9) 


From the geometry of the swirler At is 0.366 in^. Thereby, the needed diameter for the swirling flow 
case is 0.418 in. The actual jet being used has an inner diameter of 0.43 in. with an outer diameter of 03 
in. For Configuration 10 the calculated diameter is 0.591 in., with actual inner and outer diameters of 
0.591 in. and 0.75 in., respectively. The primary jet is shown in Figure 3.3.2-4. 

To determine the losses through the rig, a single swirler rig was built. This rig is shown in Figures 33.2- 
5 and 33.2-6. The test section was 3.0 in. square and extended 10 duct lengths downstream. Pressure taps 
were installed on the swirler annular pipe, primary jet, and main duct downstream of the mixing region. 
Velocity measurements were taken with a hot wire downstream of the mixing region before entrance 
into the plenum. From these measurements, the loss coefficients were calculated for the primary jets, 
annular pipe, and swirler. A summary of the measured losses is shown in Figure 33.2-7. In this figure K 
is the loss coefficient in reference to the mean annular or primary jet velocity. Also, the measured val- 
ues of static pressure and test section velocity are indicated. 

33.2.1 Data Acquisition 

The data acquisition system is shown in Figure 33.2-8. A signal from the PM tubes contains the 40 Mhz 
signal from the Bragg cell plus or minus the Doppler signal. PM signal is greater than 40 Mhz if the 
flow is into the direction of fringe movement and less than 40 Mhz if the flow is with the fringe move- 
ment. 


Since the signal frequency is higher than that allowable for the processors, it is sent through a fre- 
quency mixer. The mixer heterodynes this signal with a selected mix frequency ranging from 31 to 49 
Mhz. 

Output from the mixer is sent to the TSI counter processors, one for each color line. The processor ini- 
tially provides gain and high and low pass filtering to eliminate noise from being processed. Each par- 
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tide generates a burst as it passes through the probe volume, and at the onset of a burst a timer is 
started. If the burst contains a set minimum number of cycles, the clock is stopped. If the burst contains 
less than the required number of cycles, the burst is considered invalid and the processor resets to accept 
another burst. This step eliminates particles that pass along the edge of the probe volume which pro- 
vide a signal burst with a small number of cycles. For this experiment the cycles/bursts were set at 1 /2. 
Finally, the time for the passage of the cycles/bursts is compared to that for 1/2 cycles/bursts. If the 
two values do not compare within a certain percentage the sample is discarded. A 7 % comparison was 
used for the data acquisition. This step eliminates signals that would be generated when two particles 
cross through the probe volume at the same time. 

The data are then sent to the DEC 11/23 minicomputer. Besides taking the data, the computer also con- 
trols the probe volume movement. The laser table is movable in two directions via computer controlled 
stepping motors. The third direction of movement is accomplished by movement of the field lens, and is 
also done with a stepping motor. This provides full 3-D movement of the probe volume in the flow 
field. The axis movement sensitivities are listed below. 

• horizontal table movement: 1/1000 in. per pulse 

• vertical table movement: 1/22574 in. per pulse 

• field lens movement: 1/1000 in. per pulse 

To accommodate the use of forward scatter, as mentioned previously, a sliding boom was used between 
the transmitting and receiving optics. The boom was bolted to the bottom of the laser table and ex- 
tended under the rig to the opposite side (Figure 3.3.2-9). Bolted to the top of the boom was a sliding 
beam with one end attached to the filed lens and the other to the receiving optics. This provided full 
3-D movement of the receiving optics with the probe volume. 

3. 3.2.2 Calibration of Annular Pipe and Primary let 

To calibrate the primary jet and swirler several setups were used. The first was a setup to determine 
the flow profile (Figure 33.2-10). A high pressure air supply and plenum were used to blow air through 
one of the primary jets. At the exit, a total pressure probe was scanned across. A stepper motor driven 
traverse was used to move the probe, and a dial gage was used to measure position. From this, the pipe 
flow profile was determined. 

The measured profile is shown in Figure 3.3.2-11. This profile fits well the one-sixth turbulent profile 
presented by Schlichting (1979). Thus for the rest of the measurements a one-sixth turbulent profile is 
assumed. 

The setup used for the mass flow calibration of the primary jet is shown in Figure 3.3.2-12. This is the 
same rig used for the LDV measurements. The only difference is that the headplate for the swirler is 
replaced with one to hold the primary jet. Also, a total pressure probe is mounted on the jet centerline 
over the second pressure tap. Thus, the jet centerline pitot pressure is measured. Figure 3.3.2-13 is the 
plot of centerline pitot pressure against the pressure drop across the primary jet pressure taps. From 
Bernoulli's equation: 



( 10 ) 


solving for V 



(11) 
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Assuming a one-sixth profile the mean jet velocity is 0.791% of the jet centerline velocity (Schlichting, 
1979). 


V = 0.791, 


2AP 


Pitot 


rearranging Equation 12 


( 12 ) 


(0.791) 2 AP Pitot =jpV 2 (13) 

The relationship between the mean jet velocity and the static pressure drop across taps 1 and 2 is given 
by: 


AP 12 = K^pV 2 (14) 

Where K is the loss coefficient between the pressure taps. Using Equations 13 and 14: 

(0.791) 2 AP Pilot =^pV 2 =^ (15) 


solving for K 


AP12 

(0.791) 2 AP Pilot 


(16) 


APi2/APpitot * s the inverse slope of the plot in Figure 33.2-13. A least squares routine was used to cal- 
culate the slope. For this plot, the slope was calculated to be 3.113 and a loss coefficient of 0.513. 

With this setup it is easy to measure the losses in the other portions of the primary jet by switching 
pressure tap hook ups. Figures 3.3.2-14, 3.3.2-15, and 3.3.2-16 are sketches of the other pressure trans- 
ducer configurations. The corresponding data plots are provided in Figures 3.33-17, 3.3.2-18, and 3.3.2- 

19. 


A similar procedure was also adopted for the annular pipe. Figure 3.3.2-20 is a diagram of the setup 
used for the annular pipe mass flow calibration. Figure 3.32-21 is a plot of the measured data. The re- 
maining data plots for the annulus are provided in Figures 33.2-22, 3.3.2-23, and 3.3.2-24. These plots 
are for the same setups as Figures 3.3.2-17, 3.3.2-18, and 3.3.2-19, except that these are for the annular 
pipe instead of the primary jet. 

Finally, the same tests were run on the annular pipe with the swirler installed. However, since the 
pitot probe could not be inserted through the swirler vanes, the required pressure drops were measured 
against annular pipe pressure drop. That is, the pressure difference between taps 1 and 2 were mea- 
sured. This is easily referenced to centerline pitot pressure through Figure 3.3.2-21, which is a plot of 
centerline pitot pressure versus annular pipe pressure drop. The pressure transducer arrangement for the 
measurement of swirler loss is shown in Figure 3.3.2-25. This data is shown in Figures 33.2-26, 33.2-27, 
and 3.3.2-28. 
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3.3.2.3 Inlet Conditions 


To establish the swirler inlet flow, a finely detailed scan was taken in front of the swirler at 0.050 in. 
downstream of the exit The spatial grid used is shown in Figure 3.3.2-29. A 39 x 20 point grid was 
taken with 0.015 in. in between each nodal point, for a total of 780 spatial points. For each spatial 
point, the mean and rms were computed from the 3000 samples taken. In all, a total of 9.36 million data 
samples were taken of the swirler inlet. 

Referring to Figure 3.3.2-29, at each spatial point the laser was able to measure two components of ve- 
locity. The two components are the x component (into the plane of the paper) and y component (towards 
the top of the page). In order to get all three components, two scans were taken and the use of symmetry 
imposed. Scan 1 was able to measure u and v velocity and scan 2 was able to measure u and w compo- 
nents. Thus all three components were measured, with the u velocity being measured twice. 

Figure 3.3.2-30 contains cascade plots of the x component of velocity for scans 1 and 2. In general, the ax- 
ial velocity is lower at the hub and has maximum values at the outer flow radius. However, there is 
considerable discrepancy in the u velocity measured by the two scans. Referring to Figure 3.3.2-31, the 
maximum axial velocity varies by 40 ft/sec between scans 1 and 2. 

It is expected that there would be some spatial averaging between scans 1 and 2. The probe volume has 
a certain length and width, and for each spatial point the measured velocity is an average over the 
probe volume. In scan 1 the probe volume dimension is only 0.001 in. along the y axis, where for scan 2 
the probe volume dimension is 0.015 in. along the y direction. Thus, the spatial resolution in the y di- 
rection is considerably finer for scan 1 than scan 2. This would account for a slight variation between 
the two scans, but not a 40 ft/ sec difference. 

Plots of v and w velocities are provided in Figures 3-3.2-32 and 3.3.2-33. The v component converges to a 
relatively constant value near, or in front of, the hub. This is a strong indication of the back flow being 
drawn into the swirling exit flow. Also, the w component converges to near-zero in front of the hub, in- 
dicating low swirl in this region. As expected, at the left side of the plane (low z position) the v veloc- 
ity is high, while the w velocity is low, with the values changing as z increases. The combination of v 
and w velocities are plotted in vector form in Figure 3.3.2-34. Here, the effect of the swirl is clearly ev- 
ident. 

As mentioned earlier, rms values were also measured. Figure 3.3.2-35 is a diagram of the laser orienta- 
tion. The green and blue beams are orthogonal to each other at approximately 45 deg angles to the 
faceplate. The rms values have been nondimensionalized with the local mean velocity for that partic- 
ular beam. Thus, Figures 3.3.2-36 through 3.3.2-39 are plots of the turbulence intensity ratio or 
rms/ I mean I velocities. 

Figures 3.3.2-36 and 33.2-37 are plots of intensities for scan 1. The turbulence intensity is lowest in front 
of the annular hub or the outer flow radius. At the edge of the inner flow radius, the intensities get 
very large or, for the green beam, tend to infinity. Here the flow is very near zero and even slight fluc- 
tuations make the turbulence intensity very large. 

For scan 2 the results are similar (Figures 3.3.2-38 and 3.3.2-39). The turbulence intensities are lowest at 
the outer flow radius and tend to infinity at the inner flow radius. However, even at its lowest the tur- 
bulence intensity is still upward of 20 percent. 

3. 3.2.4 Downstream Data 

Mean and rms velocities were also taken at four downstream locations. These locations were at 0.5, 1.0, 
1.5, and 2.0 duct heights downstream of the swirler exit. For each downstream location a 26 x 31 point 
grid was used for a total of 806 spatial points, with a spatial resolution of 0.050 in. As with the inlet 
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conditions two scans were taken and symmetry invoked. The grid and scans taken are shown in Figures 
3.3.2-40 and 3.3.2-41. In all, a total of over 38 million data samples were taken downstream of the 
swirler exit. The laser orientation for each scan is shown in Figure 3.3.2-42. 

Figures 3.3.2-43 through 33.2-46 are plots of the velocity components for the first downstream location 
(x = 1.5 in). Referring to Figure 3.3.2-43, the u velocity is greatest along the walls and in the comers. 
That is, in the center of the flow there is a strong back flow, with the bulk downstream flow along the 
walls and in the comers. Referring to Figure 33.2-44 there is very good agreement between the two dif- 
ferent values of u velocity obtained from the different scans. Looking at the centerline scan for z = 13 
in., the extent of the back flow region can easily be seen. Only the half-inch nearest the wall has a 
positive u velocity. The rest of the channel is in back flow. 

The extent of the swirling flow can be seen quite clearly in Figures 33.2-45 and 3.3.2-46, and also in Fig- 
ure 3.3.2-47. For the centerline scan of w velocity at y = 2.75 (Figure 33.2-46) the w component is twice 
as strong as the u component for the same position. Figure 3.3.2-47 is a vector plot of the v and w compo- 
nents. Here, the strong swirling flow is clearly evident. Also, there is a slight asymmetry in the flow. 
That is, the apparent center of the swirl is not at the geometric center of the channel. The recirculation 
region takes the shape of the square duct at the 1.5 in. location. 

Figures 33.2-48 through 3.3.2-52 are plots for the second downstream location at x = 3.0 in. Referring to 
Figure 33.2-48, the u velocity for scan 2 shows a strong velocity deficit in the comer. This could possible 
be due to a comer vortex. However, it is uncertain if this is actually the case. It is only possible to mea- 
sure two components of velocity at these positions. Thus, a full velocity profile for these positions is 
unavailable. 

Comparing Figure 3.3.2-39 to Figure 3.3.2-44 for the centerline scan, the point of zero crossing has moved 
in towards the center of the duct. Thus, the size of the area of back flow has decreased considerably. 
Also, the magnitude of the peak velocities has decayed. For the centerline position, the maximum ve- 
locity near the wall has dropped from 40 ft/sec to 15 ft/sec. Although the strength of the flow along 
the walls has decreased dramatically, the strength of the recirculation flow has only decreased 
slightly from -25 ft/sec to -20 ft/sec. 

Comparing Figure 3.3.2-46 to Figure 3.3.2-51, the W component for the centerline scan has decreased 
only slightly. Thus, the swirl component has only decreased slightly. 

Figures 3.3.2-53 through 3.3.2-57 are plots of velocity intensity for x = 4.5 in. or 1.5 duct heights down- 
stream from the swirler exit. Referring to Figures 3.3.2-53 and 33.2-55, again there is a sudden change 
in velocity at the comers. 

Figures 3.3.2-58 through 33.2-62 are plots for x = 6.0. The plots at x = 43 and 6.0 are very similar to 
those for x = 3.0. The only difference is that there is a slight decay in the magnitude of the velocities 
but the flow profiles remain rather constant. 

Thus, the flow goes through some radical changes in the first duct height (3.0 in.). After this, the flow 
profiles remain rather constant, with only slight decay of the magnitude of the velocities. Referring to 
Figure 3.3.2-63, these are contour plots of u velocity with areas of negative or back flow showing up as 
dotted areas. From these plots the size of the recirculation region is clearly visible. At x = 1.5, the ef- 
fect of the square duct can easily be seen in the squared-off back flow region. By x = 3.0 or 1.0 duct 
height, the area of recirculation has diminished considerably and has become circular in shape. For x = 
4.5 in. and 6.0 in., the size of the recirculation region diminishes very little, if at all. This is also seen 
in Figure 3.3.2-64, which is a vector plot at the centerline or z = 13 in. position. 

In doing some flow visualization, it was later discovered that there are two exit flow regimes. These 
are shown schematically in Figure 33.2-65. At low mass flow rates, the exit flow expands outward at 


204 


an angle of considerably less than 90 deg, and is not attached to the headplate. At high mass flow 
rates the exit flow remains attached to the headplate all the way out to the side wall. With moder- 
ate flow rates a bistable flow situation is setup. That is, the flow will oscillate between the two 
regimes. 

It is believed that this is the cause of the velocity discrepancy indicated in Figure 3.3.2-31. Also, evi- 
dence of this phenomenon can be seen in Figure 3.3.2-56. For the v velocity at z = 0.3 in. there is a very 
definite jump in the velocity profile when the exit flow switches from one regime to the other and back 

again. 

33.2.5 Turbulence Quantities 

Figures 3.3.2-66 through 3.32-69 are plots of turbulence intensity for the first downstream location. As 
with the inlet condition data, when the mean velocity goes towards zero the turbulence intensity goes 
to infinity. In areas where this is not the situation, the turbulence intensity varies from 35 to 100%. 

Figures 3.3.2-70 through 3.3.2-73 are plots for the second downstream location at x=3.0 in. Referring to 
Figure 3.3.2-70, the U velocity for scan 2 shows a strong velocity deficit in the comer. This could possi- 
bly be due to a comer vortex, however it is uncertain if this is actually the case. It is only possible to 
measure two components of velocity at these positions, thus a full velocity profile for these positions is 
unavailable. 

Figures 3.3.2-74 through 3.3.2-81 are plots of turbulence intensity for x = 4.5 in. and x = 6.0 in. The plots 
at x = 4.5 and 6.0 in. are very similar to those at x = 3.0 in, the only difference is that there is a slight 
decay in magnitude but the flow profiles are similar. 
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Figure 3. 3. 1-1. Photograph of swirling jet apparatus. 
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Figure 3.3.1 -2. Concentration fields for N=0.00: 


(a) 1 frame (b) 127 frames. 
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Figure 3.3.1-3. Concentration fields for N=0.33: (a) 1 frame (b) 127 frames. 
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Figure 3.3. 1-4. Concentration fields for N=0.66: (a) 1 frame (b) 127 frames. 
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Figure 3.3.1-5. Concentration fields for N=1.00: (a) 1 frame (b) 127 frames. 
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Figure 3.3.1-6. Concentration fields for N=1.33: (a) f frame (b) 127 frames. 
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Figure 33 . 1 - 7 . Concentration fields for N=1.66: (a) 1 frame (b) 127 frames. 
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Figure 3.3.1 -8. Concentration fields for N=2.00: (a) 1 frame (b) 127 frames. 
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Figure 3.3.1-9. Relationship between spread angle and swirl ratio. 



Figure 33.2-1. Diagram of vane swirler. 
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Figure 3.3.2-2. Drawing of annular swirler inlet. 
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Figure 3.3.2-3. Determination of primary jet diameters. 
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Figure 3.32-4. Drawing of primary jet. 
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Figure 3.3.2-5. Single swirler loss test rig. 
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Figure 33.2-6. Single swirler rig for LDV measurements. 
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Figure 33.2-7. Results of swirler loss tests. 
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Figure 3.3.2-8. Data acquisition system for LDV measurements. 
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Figure 3.3.2-9. Laser table and boom assembly. 



AIR SUPPLY 



Figure 3.3.2-10. Setup for measurement of flow profile. 
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Figure 3.3.2-11. Primary jet flow profile. 
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Figure 3.3.2-12. Pressure transducer arrangement for mass flow calibration of primary jet. 
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Figure 3.3.2-13. Mass flow calibration plot for primary jet. 
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Figure 3.3.2-14. Pressure transducer arrangement for measurement of primary jet inlet loss. 
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Figure 3-3.2-15. Pressure transducer arrangement for measurement of primary jet exit loss. 
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Figure 3.3.2-16. Pressure transducer arrangement for measurement of primary jet pressure loss. 
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Figure 3.3.2-17. Plot of primary jet inlet pressure loss. 
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Figure 3.3.2-18. Plot of primary jet exit pressure loss. 
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Figure 3.3.2-19. Plot of primary jet pressure loss. 
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Figure 3.3.2-20. Pressure transducer arrangement for mass flow calibration of annular pipe. 
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Figure 33.2-21. Mass flow calibration plot for annular pipe. 
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Figure 3.32-22. Plot of annular pipe inlet pressure loss. 
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Figure 3.3.2-23. Plot of annular pipe exit pressure loss. 
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Figure 3.3.2-24. Plot of annular pipe pressure loss. 
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Figure 3.3.2-25. Pressure transducer arrangement for measurement of swirler inlet loss. 


Static Pressure Drop Across Inlet 


0. 080 




Annular Pipe Pressure Drop Cpsi) 

TE92-2438 


Figure 33.2-26. Plot of swirler inlet pressure loss. 
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Figure 3.3.2-27. Plot of pressure loss across swirler. 
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Figure 33.2-28. Plot of swirler pressure loss. 
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Figure 3.3.2-29. Swirler inlet data grid. 
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Figure 3.3.2-30. Cascade plots of swirler inlet u velocity (scans 1 and 2). 
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Figure 3.3.2-31. Swirler inlet plots of u velocity at specific Z locations (scans 1 and 2). 
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Figure 3.3.2-32. Cascade plots of swirler inlet v and w velocities. 
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Figure 33.2-33. Swirier inlet plots of v and w velocities at specific Z locations. 
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Figure 3.3.2-34. Vector plot of swirler inlet v and w velocities. 
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Figure 3.3.2-35. Swirler inlet laser orientation. 
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Figure 3.3.2-36. Cascade plots of turbulence intensity for s wirier inlet (scan 1). 
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Figure 3 . 32 - 37 . Swirler inlet plots of turbulence intensity at specific Z locations (scan 1). 
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Figure 3.3.2-38. Cascade plots of turbulence intensity for swirler inlet (scan 2). 
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Figure 3.32-39. Swirler inlet plots of turbulence intensity at specific z locations (scan 2). 
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Figure 33.2^40. Downstream data grid (scans 1 and 2). 
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Figure 3.32-41. Combined data grid. 
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Figure 3.3.2-42. Laser orientation for downstream data. 
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Figure 3.3.2-43. Cascade plots of u velocity at x=1.5 in. (scans 1 and 2). 
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Figure 3.3.2-44. Plots of u velocity at x=1.5 in. (scans 1 and 2). 
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Figure 3.3.2-45. Cascade plots of v and w velocities at X=1.5 in. 
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Figure 3.32-46. Plots of v and w velocities at x=1.5 in. 
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Figure 33.2-47. Vector plot of v and w velocities at x=1.5 in. 
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Figure 3.3.2-48. Cascade plots of u velocity at x=3.0 in. (scans 1 and 2). 
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Figure 3.3.2-49. Plots of u velocity at x=3.0 in. (scans 1 and 2). 
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Figure 3.3.2-50. Cascade plots of v and w velocities at x=3.0 in. 
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Figure 3.3.2-51. Plots of v and w velocities at x=3.0 in. 
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Figure 33.2-52. Vector plot of v and w velocities at x=3.0 in. 
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Figure 3.3.2-53. Cascade plots of u velocity at x=4.5 in. (scans 1 and 2). 
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Figure 3.3.2-54. Plots of u velocity at x=4.5 in. (scans 1 and 2). 
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Figure 3.3.2-55. Cascade plots of v and w velocities at X=4.5 in. 
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Figure 3.32-56. Plots of v and w velocities at x=4.5 in. 
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Figure 3.3.2-57. Vector plot of v and w velocities at x=4.5 in. 
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Figure 3.3.2-58. Cascade plots of u velocity at x=6.0 in. (scans 1 and 2). 
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Figure 3.32-59. Plots of u velocity at x=6.0 in. (scans 1 and 2). 
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Figure 3.3.2-60. Cascade plots of v and w velocities at X=6.0 in. 
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Figure 3.3.2-61. Plots of v and w velocities at x=6.0 in. 
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Figure 33.2-62. Vector plot of v and w velocities at x=6.0 in. 
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Figure 3.3.2-63. Contour plots of u velocity indicating regions of reverse flow. 
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Figure 3.3.2-64. Vector plot at Z=1.5 in. 
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Figure 3.3.2-65. Sketch of swirler exit flow regimes. 
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Figure 3.3.2-66. Cascade plots of turbulence intensity at x=1.5 in. (scan 1). 
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Figure 3.3J2-67. Plots of turbulence intensity at x=1.5 in. (scan 1). 
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Figure 3.3.2-68. Cascade plots of turbulence intensity at x=1.5 in. (scan 2). 
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Figure 3.3.2-69. Plots of turbulence intensity at x=1.5 in. (scan 2). 
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Figure 3.3.2-70. Cascade plots of turbulence intensity at x=3.0 in. (scan 1). 
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Figure 3.3.2-71. Plots of turbulence intensity at x=3.0 in. (scan 1). 


294 


■i Z * I » - 

rf-r" S? !s. 

fffirrffUHto,,.... 

rrfnUis!»?i. 



s 

rr 

C\1 

OJ 

05 

LU 


a I I j r 

* mi A f4 4 

M3» / 9« 


295 


Figure 3.3.2-72. Cascade plots of turbulence intensity at x=3.0 in. (scan 2). 
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Figure 3.3.2-73. Plots of turbulence intensity at x=3.0 in. (scan 2). 
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Figure 3.3.2-74. Cascade plots of turbulence intensity at x=4.5 in. (scan 1). 
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Figure 3.3.2-75. Plots of turbulence intensity at x=4.5 in. (scan 1). 
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Figure 3.3.2-76. Cascade plots of turbulence intensity at x=4.5 in. (scan 2). 
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Figure 3.3.2-77. Plots of turbulence intensity at x=4.5 in. (scan 2). 
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Figure 3.3.2-78. Cascade plots of turbulence intensity at x=6.0 in. (scan 1). 
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Figure 3.3.2-79. Plots of turbulence intensity at x=6.0 in. (scan 1). 
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Figure 3.3.2-80. Cascade plots of turbulence intensity at x-6.0 in. (scan 2). 
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Figure 3.3.2-81. Plots of turbulence intensity at x=6.0 in. (scan 2). 
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3.4 FINAL RIG - INSTRUMENTATION AND DATA REDUCTION 


Using the experience gained with the initial rigs, the swirling jet and the single swirler rig, a new air 
test rig was designed and constructed. The details of the final apparatus are described and the instru- 
mentation, data analysis, and error analysis are presented in the next subsections. 

3.4.1 Final Experimental Apparatus 

Figures 3.4.1-1, 3.4.1-2, and 3.4.1-3 are three-view drawings of the three configurations of the combustor 
model. The annular combustor model was designed to simulate the cold flow characteristics of current 
gas turbine combustors. Five cells were used in its construction in order to allow for any interaction be- 
tween cells and to remove wall effects. The primary zone of a combustor usually has three flows enter- 
ing this region. These flows consist of air from the discharge of a compressor through swirlers in a dome 
wall, fuel injection through the center of the swirler, and primary air jets normal to the mean flow in 
the combustor. No fuel injection was used in this research. 

The cross section was 15 in. wide and 3 in. high. The test section extended downstream a distance of 30 
in. Each cell contained one pair of primary jets, located 13 in. downstream of the endplate, mounted 
normal to the flow from the annular jets. 

Figure 3.4.1-4 shows a drawing of both annular and primary jets. Both the annular and primary jets 
were constructed of aluminum, and were designed with a contraction in order to provide a flat velocity 
profile upon entrance into the test section. The annular jet had an inlet diameter of 2343 in. and diame- 
ter of 1.459 in. at the test section entrance. The center portion of the annular jet had a diameter of 1.084 
in. at the test section entrance. This results in an annular gap of 0.1875 in. and entrance area of 0.7489 
in^. The primary jets were 1.375 in. in length with an inlet diameter of 1.0 in. and diameter of 0.43 in. at 
the test section entrance. This results in an inlet area of 0.1452 in^. for each primary jet. 

In order to damp out disturbances and allow uniform mixing of the seed particles before entrance into the 
test section, tubing was attached to the primary and annular jets. Tubes extending 10.5 in. above the ma- 
chined jets were connected to all 10 primary jets. Preliminary flow visualization indicated that these 
tubes helped reduce unsteady displacement of the jet stagnation point. Tubes extending 9.75 in. in front 
of the annular jets were used mainly to provide adequate length to allow mixing of the smoke. Figure 
3.4.1-5 shows a drawing of the model when annular and primary jets are attached to the model. 

The exit of the combustor model was connected to a plenum chamber, which was connected to a fan. A 
centrifugal blower was used to draw air through the test section. Previous work preformed on this 
model indicated that large-scale feedback from the fan existed. To eliminate these large disturbances 
screens and foam were placed in the last 93 in. of the combustor model. 

3.4.2 Instrumentation 


3.4.2. 1 Laser Doppler Velocimeter 

The LDV system used in this experiment incorporated a Spectra Physics model 165 Argon ion laser. It 
was configured to measure a single component of velocity using the wavelength of 4880 Angstroms, blue. 
Figures 3.4.2-1 and 3.4.2-2 show schematics of the LDV system used. 

The transmitting optics used were made up of the DISA 55X components. In the first stage of the trans- 
mitting optics, the beams entered a beam waist which was used to adjust the beam diameters to a mini- 
mum at the probe volume. This was necessary because when beams have a minimum diameter the 
wavefronts are planar. The interference pattern produced by the intersection of planar waves produced 
an evenly spaced interference pattern across the probe volume. If the wave fronts were spherical the 
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resulting fringe surfaces would be distorted and an evenly spaced interference pattern would not be pro- 
duced. Particles with the same velocity, but crossing different parts of the probe volume, would produce 
different Doppler frequencies when spherical wave fronts intersect. 

A beam splitter was used next to split the single beam into two equal intensity beams. One beam stayed 
along the optical axis, while the second beam was translated approximately 30 mm from the optical 
axis. 


The frequency of the beam along the optical axis was acousto-optically shifted by 40 MHz using a 
Bragg cell. The Bragg cell caused the fringes in the probe volume to move at 40 MHz. With fringe 
movement in the probe volume, the velocity directions of the particles crossing the probe volume could 
be detected. A frequency of 40 MHz plus the Doppler frequency resulted when particles were moving in 
the opposite direction of fringe movement, and a frequency of 40 MHz minus the Doppler frequency re- 
sulted for particles moving with the fringe movement. In addition to detection of velocity directions, 
the Bragg cell helped correct for angle or fringe bias in the velocity measurements. 

The final set of optics in the transmitting optics package was the beam translator. The beam translator 
moved the off axis beam back towards the optical axis. This was necessary in order to allow entrance 
into the beam expander. Also, adjustments in the beam translator allowed the beams to be focused to a 
common intersection point. 


Upon exit from the beam translator, a beam expander was used to increase the beam diameters by a ra- 
tio of 3.75 to 1. The purpose of increasing the beam diameters was to reduce the size of the probe vol- 
ume. With a smaller probe volume, gradient bias was reduced and greater resolution of the flowfield 
could be obtained by using smaller grid spacing. The radius and length of the probe volume are given by 
Equations 17 and 18. 



2b 0 



(18) 


where: 

b = radius of input beam 
f = focal length of field lens 
X = wavelength of laser beam 
0b = angle between beams 


A field lens, made up of two lenses, focused the two incoming beams to a common point. Focal lengths for 
each lens were i\ = 30 in. and f2=48 in. From lens theory, the effective focal length of the two lenses can 
be expressed by Equation 19. 


f = 



(19) 


Therefore, the field lens had a focal length of 18.46 in. Table 3.4-1 lists parameter settings of the LDV 
system used for this research. 
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Table 3.4-1. 

LDV parameter settings. 


Optical system setting s 
wavelength 
focal length of lens 
diameter of input beam 
angle between beams (average) 
length of probe volume (average) 
diameter of probe volume (average) 
frenquency shift 

Data collection setting s 
mode of operation 
filter settings 

number of fringes per signal 

comparator 

sample size 

data rate 

sample rate 

Flow system parameters 
velocity range 
seed particles 
seed particle size (average) 


488 nm = 1.601 x 10"^ ft 
18.46 in. 

0.175 in. 

5 deg 
0.056 in. 

0.0027 in. 

40 MHz 

N burst mode 
20 MHz low pass 
2 MHz high pass 
16 
1 % 

10,000 samples 

10,000 to 90,000 samples/sec 

1000 samples/sec 

-35 ft/sec S U < 35 ft/sec 
propylene glycol 
2 microns 


The receiving optics used consisted of different elements made by DISA. The optics were used in for- 
ward scatter, which ensures a stronger signal than in back scatter. Scattered light from the probe vol- 
ume was focused to a point by a receiving lens. Next, a beam expander working in reverse reduced the 
diameter of the light from the receiving lens, allowing entrance into a spatial filter. 

In order to reduce extra light from reflections, room light, or scattered light not from the probe volume, 
a spatial filter was used. This filter was the same size as the probe volume. The light passing through 
the spatial filter entered a photomultiplier. The photomultiplier was used to convert light intensities 
into electrical signals proportional to the intensity of the light. The output of the photomultiplier was 
then sent to a signal processing unit. 

3. 4.2.2 Data Acquisition System 

The output signal from the photomultiplier contains a signal that has been shifted by 40 MHz plus or 
minus the Doppler frequency from a measurement. This signal must first be frequency shifted to allow 
the processing electronics to handle the signal. A DISA 55N12 frequency mixer was used to mix the sig- 
nal with a frequency between 31 MHz to 49 MHz. See Figure 3.4.2-3 for a schematic of the data acquisi- 
tion system. 

After frequency mixing, the signal was sent to a TSI model 1980 counter-type signal processor. This sig- 
nal first entered a model 1984 input conditioner where high and low pass filtering occurred to remove 
the pedestal or DC component from the signal and to remove any high frequency noise. In addition, gain 
adjustment and minimum cycles/burst settings were set in the input conditioner. The cycles/burst setting 
was used to determine the number of fringes a particle must traverse in the probe volume before that 
signal is accepted by the processing electronics as a valid data point. A setting of 16 cycles/burst was 
used for all measurements. This eliminated data from particles that passed through the edge of the 
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probe volume and crossed only a few fringes. The filtered output jack was connected to a Hewlett 
Packard spectrum analyzer to monitor the Doppler frequency of the signal. 

Output of the input conditioner was sent to a model 1985 timer for further analysis. A percent compari- 
son setting is used to make a comparison between the time for N cycles and N/2 cycles. Measurements 
were made with a 1% setting, which ensured an actual 0.8% comparison between N and N/2 cycles (the 
time for 16 and 8 cycles must be within 0.8% for the data point to be accepted). Additionally, settings 
for the exponent and setting manual/ auto ranging on the exponent was available. The data ready jack 
provided a 0.6 microsecond pulse every time a new data point is validated. This was connected to a 
counter to determine the data rate. A monitor output jack is connected to a Tektronix 215A oscilloscope 
in order to observe validated data. 

Data were output to a PC's Limited microcomputer through digital outputs on the timer unit In addi- 
tion to processing the data, the computer controlled the position of the probe volume within the test rig 
via a Metrobyte Das20 A/D board. Table movement commands were carried out by three stepper motors 
allowing 3-D movement of the probe volume. A laser table contained two directions of movement, 
where the third direction of movement was accomplished by movement of the field lens. Probe volume 
position in the test rig was read into the computer by using an Acu-Rite III digital readout system. A 
scale attached to each axis of the table provided position information. Each scale had an accuracy of 
0.0005 in., which facilitated accurate repositioning of the probe volume for repeated runs. 

Velocity bias was caused by the arrival rate of particles in the probe volume being dependent on the 
flow velocity. A burst mode signal processor was used in this research, which operated on the signals 
generated by single particles passing through the probe volume. Therefore, particles passed through 
the probe volume more often when the flow velocity was high and less often when the flow velocity 
was low. The result of simply averaging these signals would have caused the mean velocity to be bi- 
ased high and the rms velocity to be biased low. 

Many different methods are presently available to correct for velocity bias. Equal time sampling was 
used to eliminate the velocity bias for this research. Roesler (1980) presents a study done on velocity 
bias and shows that this velocity bias may be eliminated by using equal time sampling. Data rates of 
100 times the sampling rate were recommended to ensure equal time sampling, however it was shown 
that data rates of 10 times the sampling rate provided equal time sampling. Data rates between 10,000 
to 90,000 samples per second were present for this study, with average rates of 45,000 samples per sec- 
ond. A sampling rate of 1,000 Hz was used providing data rates 10 to 90 times the sampling rate. With 
these data rates, equal time sampling should have existed between samples. A total of 10,000 samples 
were used to determine mean velocities for each point. 

3. 4. 2.3 Velocity Measurements and Statistics Calculation 

Since only a one component LDV system was used to determine U, V, U', V’, ITV’, and K, measurements 
in the coordinate directions, only to determine these quantities would be in error. Measurements needed 
to be taken at +/- 45 deg to the coordinate axes in the plane of each component. By taking measurements 
at +/- 45 deg both the U' and V' velocity fluctuations contributed to the instantaneous Doppler fre- 


quency. 


Referring to Figure 3.4.2-1 velocity measurements at +/- 45 deg were 




( 20 ) 



( 21 ) 
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adding Equations 20 and 21, gives 



Subtracting Equation 21 from 20 

U^-u^v 
V = ^(u us -U_ t5 ) 


The instantaneous velocity, U, is made up of a mean term, U, and a fluctuation term, U'. 
erage of the square of Equation 20 is 


+ M =I((U + U')+(V + V')) 
expanding Equation 25 results in 


U*45 + U+45 2 + - — 


(u + U')(U + U') + ( V + V')(V + V')+ 2(u + U')(v • 


||u 2 + U' 2 + 2UU' + v 2 + V' 2 + 2W' + |uV + Uh/^ UV' + VI?) j 


Since 

U ^ 2= ^( U + V ) U ^5 = ^( U ' + V ') 

The average of the product of these terms is 

2U^UUs = UU' + UV' + WF + W’ 

Therefore, terms on the right side of Equation 27 cancel resulting in 
UmJ 2 + UUs 2 = j(u 2 + U 72 + v 2 + V 72 + 2UV + 21TV 7 ) 
Squaring the mean of Equation 20 results in 

LU 2 =l(u 2 + 2 Uv + v 2 ) 

Subtracting Equation 31 from Equation 30 results in 


(22) 

(23) 

(24) 

(25) 

The time av- 

(26) 

~i- 

J (27) 

(28) 

(29) 

(30) 

(31) 
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(32) 


U^ 5 2 =-i(u ,2 + 2U'V' + V' 2 ) 

A similar result can be achieved by using the -45 deg component. 

U^ = I(u^-2U 7 V 7 +V' 2 ) (33) 

By subtracting Equation 33 from Equation 32 a relation for the Reynolds shear stress term in the xy plane 
results. 


(34) 

By adding Equation 33 to Equation 34 a relation for the rms velocities in the two coordinate directions 
and the +/- 45 deg directions results. 


U' 2 + V' 2 = U^s 2 + Ul4 S 2 (35) 

Since measurements have been made along +/- 45 deg directions and the horizontal coordinate direc- 
tion, U, an expression for the rms velocity in the vertical coordinate direction can be found. 


V' 2 = U^s 2 + UV.5 2 - U' 2 (36) 

The 2-D turbulent kinetic energy term can now be calculated by 

K = (J’ 2 + V’ 2 (37) 

The results of Equations 25, 34, 36, and 37 are used to calculate the quantities of interest. 

In order to allow measurements at +/ - 45 deg, the transmitting optics had to be rotated about the optical 
axis. This rotation limited the portion of the rig that could be measured due to beam interference with 
the edges of the rig resulting in no probe volume formation. Vertical travel of the probe volume was 
limited between 0.5 in. < y < 2.5 in. 

Orientation of the beams within the rig is seen in Figure 3.4.2-5. The beam located on the optical axis 
was aligned parallel to the lower wall of the test rig for each measurement. This orientation made the 
measurements taken not exactly along the axis intended. Figure 3.4.2-6 shows a sketch of measurements 
taken. Since the beam angle was approximately 5 deg, the cosine of the half angle, 2.5 deg, is 0.999, in- 
dicating that negligible amount of error occurs in the measurements. 

From the Doppler frequency received from the signal processor computations were done to calculate a 
mean and rms velocity. For each data point reported an average of 10,000 samples taken at 1,000 Hz 
sampling frequency was used to determine the statistics. The mean velocity was calculated by using the 
equation 


U = 


IU, 


1=1 


N 


(38) 


where 
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N = the number of samples taken 
Ui = the instantaneous velocity 

The rms velocity was calculated using 

1 

2 

(39) 


U TOS = 


N , N N 

Suf-lUtSu, 

i =1 i »1 1=1 

N 


3.4.2.4 Error Analysis 

The use of an LDV for velocity measurements has an advantage of being nonintrusive to the flowfield. 
One disadvantage is that care must be taken processing data in order to avoid statistical bias problems 
that exist with the use of an LDV system. This subsection will attempt to identify some possible causes 
of errors that accompany the use of an LDV system and estimate the uncertainty of measurements re- 
ported in this thesis. 

3.4.2.4.1 Sources of Statistical Bias 

Velocity bias occurs due to the arrival rate of particles in the probe volume being dependent on the flow 
velocity. A burst mode signal processor was used in this research, which operates on the signals gener- 
ated by single particles passing through the probe volume. Therefore, particles passed through the 
probe volume more often when the flow velocity was high and less when the flow velocity was low. 

The result of simply averaging these signals would have caused the mean velocity to be biased towards 
high values and the rms velocity to be biased toward low values. 

Many different methods were available to correct for velocity bias. This research used equal time sam- 
pling to eliminate the velocity bias. Roesler (1980) presents a study done on velocity bias and shows 
that this velocity bias may be eliminated by using equal time sampling. Data rates for the present 
study varied from 10,000 to 90,000 samples per second and a sampling rate of 1,000 Hz was used. Data 
rates remained 10 to 90 times the sampling rate. With these data rates, equal time should have ex- 
isted between samples. 

Filter bias occurs because of the tendency of real systems to have measurement efficiency dependent on 
the flow velocity. If flow velocities produce a frequency outside of the processor's frequency range the 
particle will not be processed. A biased average exists if Doppler frequencies are outside the filter set- 
tings, since these frequencies are not be processed. For this research, velocity ranges were between +/- 35 
ft/sec. The frequencies produced by these particles were well within the frequency response of the pro- 
cessor, and the filter settings were set to have a flat response in these frequencies. 

Angle or fringe bias is caused by the inability of processors to measure all velocities at all angles. 
Counter processors use a preset number of fringes that must be crossed by a particle before a signal is ac- 
cepted as valid. Therefore, some particle trajectories may not cross the required number of fringes. The 
flow angle range over which particles will be measured is called the acceptance angle. To correct for 
this bias it is desired to have the measurement system capable of measuring all flow angles encoun- 
tered. According to Edwards (1987), by using a Bragg cell with fringe velocity twice the maximum 
Doppler frequency this error can be corrected. A maximum Doppler frequency of 2 MHz was encountered 
in the research and a Bragg cell shift of 40 MHz was used. No angle bias should exist. 
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Gradient bias occurs due to having a finite probe volume size and the possible existence of several dif- 
ferent velocities in the probe volume at one time. This error was ignored in this analysis, but with an 
average probe volume length of 0.06 in. this error is very small compared to others. 

3.4. 2.5 Uncertainty Errors 

Uncertainty of measured values occurs due to finite resolution of measuring devices used to determine 
constants and settings on equipment used in the data reduction process. Equation 40 shows the relation- 
ship between the Doppler frequency, beam angle, and wavelength of laser light used to determine the 
flow velocity. 



Uncertainties existed in the calculated mean velocity due to beam angle calculation and resolution of 
the processor to provide the Doppler frequency. In addition, a sampling error existed due to a finite 
number of samples taken and an error due to rotation of the optics at +/- 45 degrees. All contributed to 
the total uncertainty of the velocities reported. These individual uncertainties were added together to 
provide a total uncertainty in reported values by the sum of squares technique, Equation 41. 


AU tot =^(AU 9 ) 2 + (AU fD ) +-(AU samp ) + (AU ral f 
3.4.2.5.1 Beam Angle Uncertainty 

The uncertainty due to the measurement of the beam angle is calculated by using Equation 42. 
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For a beam half angle of 2.5 deg and A©=0.1 deg (0.001745 radians), an uncertainty was calculated for a 
representative case at x=2 in., z=7.5 in. for the annular and primary jets case. Figure 3.4.2-7 shows the 
results of this calculation on the mean velocity. 

3.4.25.2 Frequency Uncertainty 


The uncertainty due to the processor frequency resolution was determined by Equation 43. 

Jt 


3U 

3f 


D 2 sin 




(43) 


The uncertainty in Doppler frequency comes from taking the partial derivative of Equation 44 with re- 
spect to the digital mantissa, D m . 


^ D 


NxlQ 9 
D m x 2 (n_z) 


(44) 
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where 


N = number of cycles per burst 
Dm = digital mantissa 
n = exponent 

resulting in 


3f n _ NxlO 9 
3D m Dj~2 F 2 * 


AD„ 


( 45 ) 


The resolution of the digital mantissa is 1 bit in 4096, ADm=l- Substituting the result of Equation 45 
into Equation 43 and rearranging results in 





f n 2 2 (n ~ 2) 

NxlO 9 


AD m 


(46) 


For a beam angle of 25 deg, the uncertainties due only to the frequency resolution of the mean velocity 
are seen in Figure 3.4.2-8. 


3.4.2.53 Sampling Uncertainty 

Since only a finite number of samples were taken to determine the mean velocity, a sampling error will 
exist. If an infinite number of samples were taken, the true mean would have been reported and, there- 
fore, this error would no longer exist. The 95% confidence interval for the mean velocity is 


U - 1.96 



< U T £ U + 1.96 



where 

Urms = root mean square velocity 

U = sampled mean velocity 

N • = number of samples 

U T = true mean velocity 


For the rms velocity, the 95% interval is 
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1.96 
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(47) 


(48) 


where 


U m , ST = the true root-mean-square velocity 

To determine the mean values, 10,000 samples were taken. Figure 3.4.2-9 represents the uncertainty due 
to this sampling interval for the mean and rms velocities respectively. 
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3.4.2.S.4 Rotational Uncertainties 

In the previous subsection, relations were developed relating velocities measured at +/- 45 deg to the 
coordinate axis and along one coordinate axis to calculate values for V, V’, U’V’, and K'. Uncertainties 
due to accurately rotating the beam to +/- 45 deg and 0 deg within the rig will arise in uncertainties in 
the calculated values (Durst, Melling, and Whitelaw, 1976). 

Uj = U 1 cos(a 1 ) + U 2 cos(a 2 ) (49) 

U( 2 = U 1 ' 2 cos 2 (a 1 )+ U 2 2 cos 2 (a 2 ) + 2U 3 U 2 cos(a 1 )cos(a 2 ) (50) 

Equation 49 is used to relate the mean velocities measured along lines at angles of +/ -45 deg and 0 deg to 
the axial coordinate (x axis). Equation 50 is the relation between the fluctuating velocity components of 
the velocities measured at angles to the coordinate axis. 

Equations 49 and 50 can be put into matrix form and solved for the calculated quantities V, V', and 
U'V'. 
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cos 2 (a ix ) sin 2 (a ix ) cos(a lx ) sin(a ix ) 

' U' 2 ' 

u^ 2 

= 

cos 2 (a 2X ) sin 2 (a 2X ) cos(a 2X ) sin(a 2X ) 

V' 2 

U' - 2 

u -45 


_cos 2 (a 3x )sin 2 (a 3X )cos(a 3X )sin(a 3x )_ 

2U'V' 


(51) 


(52) 


Where ai x is the angle of measurement 1 with the x axis. By making the assumption that the mea- 
surements along the x axis is correct and no correction is needed, then Equation 51 becomes 


U.45 


L U -*J 



cos(a 2X )sin(a 2X ) 

cos(a 3X )sin(a 3X ) 


U 

V 


(53) 


To determine uncertainties in V, V’, and U’V’, 0 deg is substituted into angle a-[ x and +/- 45 deg +/-0.5 
deg is substituted into angles a2 x and a3 x . These matrices can then be manipulated into equations re- 
lating calculated quantities to measured quantities. Relations in the following form result 


V = C, U+45 + C 2 U_45 

V' = C 3 U 0 +C 4 U +45 +C 5 U^5 (54) 

ur = c 6 u 0 + c 7 u^ + c 8 u^ 

These relations were then subtracted from the equations derived previously, resulting in the uncer- 
tainty due to rotating the optics. Figures 3.4.2-10, 3.4.2-11, and 3.4.2-12 show the uncertainty due to ro- 
tating the optics for V, V’, and U’V’ respectively. 

3.4.2.6 Total Uncertainty 

Total uncertainty due to all previous sources are then calculated by using Equation 41. Figures 3.4.2-13 
through 3.4.2-18 represent the total uncertainty for U, U’, V, V’, U’V’, and K. 
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3.4.2.7 Marker Nephelometrv 


Marker nephelometry is a technique by which concentration fields can be studied by seeding the fluid 
stream with passive marker particles and detecting the marker particle concentration by a light scat- 
tering, optical method. Rosenweig (1966) first developed the theory and basic method and, along with 
Hottel and Williams (1961), was the first to apply the technique. Becker (1977) gives a very thorough 
and complete analysis of the technique. 

The measured concentration field is relative in nature, not absolute. The actual concentration at a point 
in the flowfield cannot be directly measured, only the concentration at that point relative to the con- 
centration at another point in the flowfield. This method requires the use of some type of light source 
(laser, lamp, etc) to illuminate the flowfield of interest, particles to scatter the incident light, and a 
detection device (photomultiplier tube, SIT, etc) to collect the scattered light. In some flowfields, nat- 
urally occurring particles scatter adequate light for the method to be successful but in most cases some 
type of light scattering seed particle must be added to the flow. 

Ideally, the light scattered by marker particles is a linear function of the number of marker particles 
residing in the measurement volume. This linear behavior depends on two criteria, independent scat- 
tering and monodisperse marker particles. The first condition is rather easily satisfied as, from Becker 
(1977), independent scattering occurs for marker particle density sufficiently small such that the center- 
to-center distance between particles is greater than 3 radii. This corresponds to a marker particle vol- 
ume fraction of 30%. For this experiment, this condition was easily satisfied. 

The second condition, the need for monodisperse particles, was not as easily met. Monodisperse parti- 
cles are all of the same size. Since nearly all particle generators created polydisperse particles, the in- 
tensity of the scattered light will not be linearly related to the number of particles in the measurement 
volume but will be a weighted function of the various particle sizes existing there. 

From Schaughnessy and Morton (1977) the logical assumption to be made is that the statistical distri- 
bution of particle sizes is the same at every point in the flowfield. This assumption will fail when the 
sample of marker particles is so small (very low concentration) as to not characterize the correct distri- 
bution of sizes. This did not create a problem in this experiment since the essential characteristics ob- 
served were dictated by the higher levels of concentration. 

Both smoke, in the air experiments, and a fluorescent dye, in the water experiments, were used in the 
course of the investigation. As a test case, a fluorescent dye (Rhodamine 6-G) was injected into the jet 
flow far upstream (22 jet diameters) of the jet exit to allow thorough mixing. The beam of a Spectra- 
Physics argon ion laser was passed through a cylindrical lens, producing a thin sheet of light aligned 
with the jet axis. The sheet of laser light was spread much wider than the measurement region to re- 
duce the effect of the Gaussian intensity variations across the laser sheet. The laser light caused the 
dye to fluoresce, such that the fluorescent intensity was proportional to the laser power and the local 
concentration of dye, for laser intensities below the fluorescence saturation level as described by 
Walker (1987). Koochesfahani and Dimotakis (1984) developed an expression for the dye concentra- 
tion at a point, i. 

C(i) = K(i) V(l ^ - (l) (55) 

where 

C ( i ) = local dye concentration 

K(i) = attenuation factor for absorption of the laser intensity by the dye 
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V ( i ) = pixel output 


D(i) = dark response of ith pixel 

The quantity fi(i) contains the effects of nonuniformity of imaging and pixel sensitivity and laser sheet 
thickness variations. These effects were experimentally determined to be negligible by observing the 
measured fluorescent intensity variations across the test section with a uniform dye concentration (see 
Appendix). The attenuation factor, K(i), was kept very near 1.0 by injecting dye at a very low rate 
(approximately 20 ppb), thus reducing the amount of dye available to absorb the laser power. The 
choice of Rhodamine 6-G for the dye also reduced the absorption. Figure 3.4.2-19 shows the mean con- 
centration profile of Fluorescein, another fluorescent dye, in which the laser light entered the jet from 
the right side of the axis. The apparent variation in concentration across the jet is due to the absorption 
of the laser intensity by the Fluorescein. This did not occur with the Rhodamine 6-G, which will be 
demonstrated by the results later in this subsection. To determine absolute concentration levels, the sys- 
tem would have to be calibrated to provide concentration for a given local laser power. However, only 
relative values of concentration were required for this study. Therefore, the local dye concentration is 
the difference between the pixel output and the dark response. The dark response was obtained by 
recording the camera output with no dye in the test section. 

Figure 3.4.2-20 shows a schematic diagram of the test equipment setup. The fluorescent intensity was 
recorded through a Panasonic digital video camera with a pickup element consisting of a 499 x 574 ma- 
trix of charge coupled devices, which will be referred to as segments of the pickup element. The camera 
was operated in the strobe effect shutter mode, which records each field in 1/1000 sec instead of the 
normal 1/60 sec. This reduces the time averaging across each field. The fields are still 1/60 sec apart to 
remain compatible with standard video cassette recorders (VCRs). Two consecutive fields are inter- 
laced to produce a frame. Figure 3.4.2-21 illustrates how the camera acquires the fields and how they 
are interlaced. The effect of the 1/60 sec between fields can be seen by plotting the odd and even pixels 
separately on a vertical scan across the image. This is shown in Figure 3.4.2-22, where the concentra- 
tion values for a single frame are plotted. The largest changes between fields occur at the jet edges, as 
expected. This has no effect on the mean and rms calculations. 

The camera output was recorded by a GE VCR onto standard 1 /2 in. tape. The recorded images were dig- 
itized through an IBM PC/AT equipped with Data Translation’s DT-IRIS image processing software 
and hardware. The digitized image is the result of assigning a value from 0 to 255 (8 bits) to each pixel 
(480 lines by 512 pixels) based on the recorded intensity. Mean and rms values at selected pixel loca- 
tions can then be computed from the digitized images. These quantities were calculated from 50 frames, 
each spaced approximately 1.4 seconds apart. The time between processed frames is due to the com- 
puter processing time. Only 50 frames were used because extensive dye recirculation within the test sec- 
tion was encountered after approximately 70 seconds. Mean and rms values were calculated as follows. 



(56) 



(57) 


N-l 


where: 
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c 

c' 

N 

Ci 


= relative mean concentration 
= rms concentration; 




= number of frames 

= instantaneous relative concentration (V(i) - D(i)) 


Visualization images are produced by mapping the intensity levels to colors, which are combinations of 
red, green, and blue, to provide pseudocolor representations of the concentration field. Pink indicates 
the highest concentration and dark blue corresponds to the lowest. These color enhanced images allow 
visualization of the concentration distribution. Figure 3.42-23 shows the 127 frame average pseudo- 
color image of the single jet dark response. 

It is desirable to have the concentration measurement system be independent of the flow velocity, espe- 
cially for this study since no velocity measurements were made. To verify that the measurement system 
was independent of velocity, the fluorescent light energy incident on a pickup segment was calculated. 
The fluorescent light intensity received by a single segment of the pickup element is 

I(t) = [l 0 (t)*e] j C(x, y, z, t)dxdydz (58) 

p.v. 

Where several quantities are defined as 

I( t) = fluorescent intensity at the pickup segment 
I()(t) = incident laser light intensity 
g = fluorescent efficiency 
C = dye concentration in the probe volume 
p.v. = probe volume 

Let 

|C(x,y,z,t)dxdydz = C(t) 

p.v. 

Then the total fluorescent light energy incident on the pickup segment during the camera exposure 
time,t e xp is 

l "p 

Ejot = Ap* J I(t)dt = Ap* J[l 0 (t)*e]C(t)dt (59) 

o o 


where 


Ap,. = pickup segment area 

Assuming that the incident laser intensity is constant with time, the total fluorescent light energy can 
then be written as 

Etot = (io* e )^ t expA p6 (60) 

where 
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C = C time averaged over texp 

From the above expression, it is apparent that the concentration measurement system is independent of 
velocity. It is also shown that the concentration measurement is averaged over the exposure time and 
over the probe volume. 

Table 3.4-II shows the sampling time and exposure time requirements for measuring various quantities. 
There are no sampling time requirements for determination of the mean, and the only requirement for 
rms determination is that the exposure time should be less than the flow time scale. To maintain the 
influence of the smallest scale fluctuations, this would require that the exposure time be less than the 
Kolmogorov time microscale. The spectrum can be determined if the time between samples is also small 
enough. 


3-4. 2. 8 Test Case - Single Turbulent let 

The measurement system was tested on a single round turbulent jet because concentration measurements 
have been made for this configuration by others. This demonstrated the validity of the measurement. 

The experimental configuration consisted of a jet of water issuing from a 0.43 in. inside diameter tube, 
with a length-to-diameter ratio of 25. The water flow was supplied to the tube from a constant head 
tank. The tank level was maintained constant by a pump which supplied water to the tank and an 
overflow pipe which discharged to the holding tank. The jet discharged vertically into a clear plexi- 
glass water reservoir. The apparatus is shown in Figures 3.4.2-24 and 3.4.2-25. Mass flow rate was de- 
termined by weighing a sample of water collected over a timed interval. The mass-averaged jet veloc- 
ity was 1.78 ft/sec, giving a Reynolds number based on jet diameter of 5300. The dye mixture, which was 
3.8 ppm dye, was injected at the rate of 0.00058 lb/sec, giving a dye concentration upstream to the jet 
exit of 20 ppb. 

For this test case, the receiving optics were arranged such that the length of the probe volume was 0.006 
in. This was determined by counting the pixels in the digitized image corresponding to a known dimen- 
sion, the jet diameter. As stated previously, the camera was operated in strobe effect shutter mode, 
which gave an effective exposure time of 0.001 seconds. 

The centerline turbulence intensity, u/U, increased through the measurement section, which extended 
5.5 jet diameters downstream. From Hinze (1975), the maximum value within the measurement region 
would occur at 5.5 diameters and be approximately 15%. Borrego and Olivari (1980) show that at 55 
diameters downstream the centerline velocity would still be equal to the exit velocity. Therefore, the 
mean velocity U is the exit velocity of 1.78 ft/sec. From the above equations, the Kolmogorov length 
and time microscales are, respectively, 0.003 in. and 0.0045 sec. This indicates that all but the smallest 
scale fluctuations are being recorded. As indicated in the previous section, this has no effect on the de- 
termination of the mean. The effect on the rms value will be minimal, since only the smallest scale 
fluctuations are being averaged across the exposure time. Since these scales are smaller than the mea- 
surement system can resolve, it is desirable to determine how much of the fluctuation energy is being 


Table 3.4-II. 

Sampling requirements for determination of various quantities. 
Quantity Sampling time required Exposure time required 

C none 

c' none 

spectrum <flow time scale 
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none 

<flow rime scale 
<flow rime scale 


recorded. Becker et al (1967), present a concentration fluctuation spectrum recorded on the jet centeriine, 
plotted against Ay(K), which is the product of the axial integral scale and the wavenumber (Figure 
3.4.2-26). They developed a linear relationship for A 

Ay = 0. 0445X (61) 


where: 

x = axial distance from the jet exit 

Based on the 1/60 sec between image fields, assume that the highest frequency which can be resolved is 
20 Hz. This would give a wavenumber of 71 ft**. The spectrum presented by Becker et al, was measured 
at 32 diameters downstream, and therefore, the current value of Ay(K) must be multiplied by 32/5.5 to 
be compatible with their spectrum. This ratioed value for Ay(K) is 3.6, giving a fluctuation energy 
level which is less than one-tenth of the energy level at low wavenumbers. Therefore, the vast 
majority of the fluctuation energy is being recorded. 

The results of this test case are presented in Figures 3.4.2-27 through 3.4.2-35, with results from other 
investigators where available. Figures 3.4.2-27 and 3.4.2-28 are visualization images of the jet concen- 
tration field. Figure 3.4.2-27 shows a single frame of data, while Figure 3.4.2-28 is the result of averag- 
ing 127 frames. 

Figure 3.4.2-29 presents the mean and fluctuation intensity c'/C data in a 3-D graphical format. The 
mean profile can be seen to flatten and spread downstream as expected. The fluctuation intensity in- 
creases outside of the peaks at the jet edges. This is due to the nearly zero rms being divided by a de- 
creasing mean value. 

Figure 3.4.2-30 presents the concentration profile near the jet exit (x/D = 0.1). The shape of the mean 
concentration profile, flat across the jet and nearly zero outside of the jet, indicates that the dye was 
thoroughly mixed with the jet flow. The fluctuation intensity has an asymmetric feature, this is due to 
the mean concentration curve. The value of C on the left side of the jet axis is slightly higher than that 
on the right side. The jet did spread faster on this side and it was visible in the flow visualization 
provided by the fluorescent dye. This was apparently due to a manufacturing defect inside the tube 
near the exit. 

Figures 3.4.2-31 and 3.4.2-32 compare results with those obtained by Rosensweig, Hottel, and Williams. 
They measured concentrations in an air jet marked with smoke particles by recording the scattered light 
from the particles. Figure 3.4.2-31 presents the mean concentration profile at x/D = 1.0. The data from 
Rosensweig et al, are asymmetric. The authors do not discuss this feature but it may be due in part to 
the fact that the projected slit width of their receiving optics was 14% of the jet diameter. This would 
make it difficult to accurately measure the concentration at the jet edges. The present data do agree 
fairly well on the left side of the jet axis. Figure 3.4.2-32 presents the fluctuation intensity profile at 
x/D = 1.5. As with the jet exit profile, there is a slight asymmetry in the profile at this location due to 
the defect in the tube. For the present data, a 95% certainty band was calculated by assuming a normal 
distribution for C, and utilizing Student's t and the % distributions. The values on the right side of the 
jet axis show good agreement with those obtained by Rosensweig et al. 

Measurements made at downstream locations of 2 and 4 jet diameters are compared with those made by 
Borrego and Olivari in Figures 3.4.2-33 and 3.4.2-34. They measured oil smoke concentrations, in a round 
turbulent air jet exhausting into ambient still air, with a light scatter meter. The radial locations are 
normalized by 5c_, the location at which the mean concentration is one-half of the centerline value. 
This allows direct comparison with Borrego and Olivari’s results, finding that the present mean concen- 
tration levels outside of the jet are higher than those measured by Borrego and Olivari. This is due to 
the recirculation of the dye within the reservoir, causing a higher background concentration. Borrego 
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and Olivari's jet exhausted into a room and therefore did not have this condition. Also included for 
comparison are concentration profiles from data recorded in air. The configuration was an air jet ex- 
hausting into still ambient air, where the light scattered by marker particles was recorded on video 
tape. The videotaped data were processed for comparison with measurements made in this study. The 
air data agree quite well with the present results, with the exception that at x/D = 2.0, the light in- 
tensity saturated the image processing system near the jet centerline. This resulted in a flat mean pro- 
file and zero rms. Most of Borrego and Olivari's rms data fall within the 95% certainty band about the 
current data. 

Figure 3.4.2-35 shows the jet half width growth along with air data and Birch et al (1978) and analyti- 
cal results from Squire and Trouncer (1944), which are taken from Forstall and Shapiro (1950). The cur- 
rent data agree well with the trend from others. As explained above, the current test hardware had a 
slight flaw, causing the jet to spread slightly faster than would be expected. 

Based on the results of this test case and the comparisons with previous investigators, the measurement 
system has been found to provide concentration measurements with satisfactory accuracy. 

3.4.2.9 Errors 

The effect of nonuniformities in the laser sheet and in the camera response were not accounted for in this 
work. In this section, those effects are investigated. Uniform concentration fields were obtained in the 
model annular combustor test section by circulating a very dilute dye solution through the test rig. The 
fluorescent field was recorded for various uniform concentration levels. Dark response data were also 
recorded with no dye in the rig. These raw data are presented in Figure 3.4.2-36. The values on the ver- 
tical axis are the levels assigned by the image digitizer. The increase at the sides of the channel are 
due to the laser sheet reflecting from the plexiglass sides. This reflection resulted in a higher laser in- 
tensity near the wall. 

Figure 3.4.2-37 presents the profiles obtained for the different concentration levels by subtracting the 
dark response from the measured intensity field. This is the method of data correction which was used 
in this study. As expected, the profiles were fairly flat across the center portion of the channel (0.5 in. 

< y < 2.0 in.). The response at the sides was apparently nonlinear since the profile drops off at each 
side. This could be due to a combination of the reflection from the plexiglass walls and nonuniformities 
in the camera response towards the edge of the image. Long, Chu and Chang (1981) suggest correcting 
for the camera's nonuniformities by recording an image of a uniform dilute concentration field and di- 
viding the measured image by this dilute response image. This is the second step of a two-step correc- 
tion which they used. The first step being the subtraction of the dark response mentioned above. 

Figure 3.4.2-38 presents the results after performing the two-step correction. The lowest concentration 
level image was used for the dilute response. The drop at the sides was still apparent and was there- 
fore probably due to the reflection from the plexiglass sides. The profiles over the center portion of the 
channel were fairly flat, but not significantly improved from those in Figure 3.4.2-37. Therefore, the 
single step correction used in this study was probably sufficient over the majority of the flow area. The 
near wall regions could possibly be improved by masking the laser sheet such that the visible sheet 
width is less than the flow channei width. The sheet must still be maintained much wider than the 
measurement area to reduce the Gaussian variation effects. 
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Figure 3.4. 1-1. Annular jets only. 
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Figure 3.4.1-2. Primary jets only. 
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Figure 3.4. 1-3. Annular and primary jets - three-view drawing. 
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Figure 3.4.1-4. Primary and annular jet - detail drawing. 
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Figure 3.4. 1-5. Head assembly - detail drawing. 
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Figure 3.4. 2-1. Receiving oplies boom support. 



Figure 3.4.2-2. LDV optics and rig setup. 
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Figure 3.4.2-3. Schematic of LDV data acquisition system. 
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Figure 3.4.2-4. Velocity measurement coordinate relations. 
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Figure 3. 4.2-5. Orientation of beams for velocity measurements. 
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Figure 3.4 2 - 6 . Actual and desired velocity component relations. 
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Figure 3. 4.2-7. U velocity error due to uncertainty in beam angle. 
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Figure 3.4.2-S. U velocity error due to frequency resolution. 
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Figure 3. 4. 2-9. U and Urms velocity error due to sampling uncertainty. 


333 


30 



y (mcnas) 


TE92-2508 


Figure 3.4.2-10. V error due to rotation of optics. 
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Figure 3.4.2-11. V velocity error due to rotation of optics. 
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Figure 3.4.2-12. ITV error due to rotation of optics. 



Figure 3.4.2-13. Total U velocity error. 
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Figure 3.42-14. Total Urms velocity error. 
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Figure 3.4.2-15. Total V velocity error. 
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Figure 3.4.2-16. Total Vrm S velocity error. 
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Figure 3.4.2-17. Total LTY’ error. 


337 


MEAN CONCENTRATION (C/C 



y (inches) 

TE92-2516 

Figure 3.42-18. Total turbulent kinetic energy error. 
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Figure 3.4.2-19. Mean concentration of fluorescent dye in round jet (x/D=0.5). 
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Figure 3.42-20. Test equipment schematic design. 


339 



0.001 sec. 


0.001 sec. 






(exposure 



FIELD #1 


FIELD #2 



(X) 


(0) 


1/60 sec. 
(sampling time) 


time 


TE92-25 20 


Figure 3.4.2-21. Temporal and spatial image acquisition scheme. 
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Figure 3.4.2-22. Instantaneous concentration profiles for consecutive image fields (x/D=1.5). 
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Figure 3.4.2-23. One hundred twenty-seven frame average dark response image of single jet setup. 
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Figure 3.4.2-24. General diagram of experimental set-up. 
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Figure 3.4.2-25. Detail of single jet test section. 
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Figure 3.4.2-26. Concentration fluctuation spectrum on jet centerline (x /D=32). 
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Figure 3.4.2-27. Single frame concentration field visualization image. 
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Figure 3.4.2-28. One hundred twenty-seven frames concentration field visualization image. 
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Figure 3.4.2-29. Three-dimensional representation of mean and rms concentration profiles. 
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Figure 3.4.2-30. Mean and rms concentration profiles at jet exit (x/D=0.1). 
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Figure 3.4.2-31. Mean concentration profile at x/D=1.0. 
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Figure 3.4.2-32. Concentration profile for rms at x/D=1.5. 
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Figure 3.4.2-33. Mean concentration profiles at x/D = 2.0 and 4.0. 
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Figure 3.4.2-34. Concentration profiles for rms at x /D=2.0 and 4.0. 
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Figure 3.4.2-35. Jet concentration half width growth. 
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Figure 3.4.2-36. Uniform concentration profiles, raw data. 


— various conc. levels 

— - DARK RESPONSE 



0.5 


1.0 


1.5 


2.0 


2 . 


m 

2 

o 

Cl 

in 

DJ 

DC 

12 

DC 

< 

Q 


>* 

in 

2 

w 

E- 

2 



y LOCATION (INCHES) 


TE92-2536 

Figure 3.4.2-37. Uniform concentration profiles, dark response subtracted. 
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Figure 3.4.2-38. Uniform concentration profiles, two step correction applied. 
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